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We generalise the Jacob-Wick helicity formalism, which applies to sequential de¬ 
cays, to effective field theories of rare decays of the type B — > Kj k (—>• Kir)£i£ 2 - 
This is achieved by reinterpreting local interaction vertices bV as a 

coherent sum of 1 —> 2 processes mediated by particles whose spin ranges between 
zero and n. We illustrate the framework by deriving the full angular distributions 
for B — > K£\l 2 and B —> K*{-^ Ktt)£i £2 for the complete dimension-six effec¬ 
tive Hamiltonian for non-equal lepton masses. Amplitudes and decay rates are 
expressed in terms of Wigner rotation matrices, leading naturally to the method 
of moments in various forms. We discuss how higher-spin operators and QED cor¬ 
rections alter the standard angular distribution used throughout the literature, po¬ 
tentially leading to differences between the method of moments and the likelihood 
fits. We propose to diagnose these effects by assessing higher angular moments. 
These could be relevant in investigating the nature of the current LHCb anoma¬ 
lies in Rk = B{B —> Kfi + [A~)/B(B —> Ke + e~) as well as angular observables in 
B -> K*n + n~. 
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1. Introduction 


Helicity amplitudes (HAs), as defined by Jacob and Wick [1], describe A —> BC (1 —> 2) 
transitions and have definite transformation properties under rotation. The key idea is that 
the angular and helicity information are equivalent to each other. Angular decay distributions 
follow (e.g. [2-4]) from evaluating the HAs with B and C in the forward direction, with the 
angular information encoded in Wigner D matrix functions, reminiscent of the Wigner-Eckart 
theorem. 

The intent of this paper is to generalise this method to decays of the type A —> ( B\B 2 )C 
which are schematically described by local interactions of the form 

H eS ~(AC) lil ... lin (B 1 B 2 y'~»'. (1) 

We do so by rewriting the 1 —» 3 decay as a sequence of 1 —> 2 processes, by inserting multiple 
complete sets of polarisation states between the Lorentz contractions of AC and B 1 B 2 above. 
This leads to a reinterpretation of the decay in terms of a sum over intermediate particles of 
spin J, where J can range from 0 up to n depending on the specific structure of the operators. 
Symbolically we may write 

n 

A(A -> (HiH 2 )C) = ^ A(A -> Bj {-+ B^C) , (2) 

j> 0 

with A denoting the amplitude. We refer to this case as the H-particle factorisation approx¬ 
imation. At the formal level, the main work is the decomposition of the Lorentz tensors into 
irreducible objects under the spatial rotation group (reminiscent of the 3+1 decomposition of 
cosmological perturbation theory for example). 

Important examples of such decays are given by the rare radiative decays B —> Kl + l~ and 
B -+ K*(—> Kir)l + £~. Besides evaluating non-perturbative matrix elements to these decays 
(e.g. [5-17]), it has become clear that it is beneficial to consider general properties of the 
amplitudes entering the angular distributions (e.g. [18-21]). Our work can be seen to be part 
of the latter category. 

We evaluate the B —> angular decay distributions within the generalised helicity 

framework developed in this paper, providing an alternative method to traditional techniques 
using Dirac trace technology [22,23]. An important consequence of the manner in which we 
derive the distribution is that it lends itself to the methods of moments (MoM), which use the 
decomposition of the distribution into orthogonal functions to obtain observables independently 
of each other. This is a complementary method to the likelihood fit to extract the dynamical 
information from the decay, and was recently studied from an experimental viewpoint in [24], 
We discuss the impact of including higher partial waves in both the (Kn)- and especially 
the dilepton-system. The latter give rise to corrections, in the form of higher moments, to the 
standard form of the angular distribution used in the literature. The sources of higher dilepton 
partial waves are higher spin operators and electroweak corrections, both of which we discuss 
qualitatively. The two sources can be distinguished by their different behaviour in higher 
partial moments. We encourage experimental investigation of higher moments from various 
viewpoints. In particular, we discuss how higher moments can be used to diagnose the size of 
QED effects in B —»• K£ + £~ (with i = e, //) and test leakage of J/T-contributions into the lower 
dilepton-spectrum. Both are of importance in view of Rk as well as the angular anomalies 
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in the low dilepton-spectrum, which have recently been reported by the LHCb collaboration 
in [25] and [26,27] respectively. 

The paper is organised as follows. In section 2 the methodology is introduced ending with 
a formal expression for the fourfold decay distribution in terms of rotation matrices and HAs. 
Specific angular distributions for B —> K^*>£ 1 ^ 2 , 1 with detailed results in appendices C (and a 
Mathematica notebook in the arXiv version) and D, are given in section 3. The method of total 
and partial moments is presented in section 4. Section 5 contains the discussion of including 
higher partial waves: a qualitative assessment of higher spin operators and QED corrections is 
presented in subsections 5.2 and 5.3 respectively. The relevance of testing for higher moments is 
emphasised in subsection 5.4. The paper ends with conclusions in section 6. Additional mate¬ 
rial, such as the leptonic HAs and a few brief remarks on A;, —> A (—> ( p , n) n) l\i 2 , is presented 
in appendices A.3 and E respectively. In appendix B we provide the kinematic conventions for 
computation of the angular distribution by the sole use of Dirac trace technology. 

2. Generalised Helicity Formalism for Effective Theories 

We first review the standard helicity formalism in section 2.1, and qualitatively apply it to 
sequential 1 —> 2 decays in section 2.1.1, specialising to the spin configuration relevant for our 
decays at the end. In section 2.2 the formalism is extended to include decays like B —> Kj k (—>■ 
Ktt)£ 1^2 described by effective field theories for b —> st\l 2 transitions. The framework can 
be straightforwardly applied to the entire zoo of semi-leptonic and rare flavour decays such 
as B s -> K*£v, B -> D^iu, D -> (ir, p)fj.fi, D —> ( 7 T,p)fji/, K —»• 7r fjp etc., and can also be 
extended to include initial particles with non-zero spin. 

2.1. The basic idea of the Helicity Formalism and its Extension 

The discussion in this section is standard and we refer the reader to [2-4] for more extensive 
reviews, as well as the pioneering paper of Jacob and Wick [1], In a 1 —> 2 (say A —> B 1 B 2 ) 
decay a particle of spin Ja and helicity Ma decays into two particles of momentum p\ and 
P 2 with helicities Ai and A 2 respectively. In the centre-of-mass frame ( p\ = —p 2 ) the system 
can be characterised by the two helicities and the direction (i.e. the solid angles 0 and (ft). By 
inserting a complete set of two-particle angular momentum states the corresponding matrix 
element can be written 

A(A ->• BiB 2 ) = (0,0, Ai,A 2 | Ja,M a ) = (ft , Ai, A 2 |j, m, Ai, A 2 ){j, m, Ai, A 2 | Ja, M a ) 

j,m 

= (6, (ft, Ai, A 2 | Ja, Ma, Ai, A 2 ) (Ja, Ma, Ai, A 2 | Ja, Ma) (3) 

'---v-''-v- 

as a product of Wigner D-functions and a HA Ai Aa . The Wigner matrix is a (2J^ + 1)- 
dimensional 50(3) representation in the helicity basis. The essence is that the distribution 
of the amplitude over the angles is then governed by the rotation matrix as a function of the 
helicities. In practice one only needs to compute the HA. 
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Throughout this work, we use non-equal leptons A £ 2 , in order to accommodate semi-leptonic decays of 
the type B —> piv as well as potential lepton flavour violation [28], motivated by the Rk measurement. 
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The process B —> J/ty(—> £ + £~)K*(—> Kir) constitutes a well-known example of a sequential 
1 —> 2 decay where the formalism can be applied [29]. The idea of this paper is to extend this 
formalism to the case where the tn^-pair emerges from a local interaction vertex Oij ~ sTibiTji 
with effective Hamiltonian H cS ~ ^2ij CijOij. This is achieved by reinterpreting the local 
interaction vertex as originating from a sum of particles whose spin depends on the number 
of Lorentz contractions between the Tjj structures. Elements of this program have appeared 
in the literature, e.g. [30] for B —> Kj £ + £~, but we are unaware of a systematic presentation 
that allows the incorporation of a generic effective Hamiltonian as well as other decay types. 


2.1.1. Helicity formalism for Bj b —> Kj k (—> K\K2)~1j 7 (—» l\ln) 

Let us consider the following sequential decay 

B Jb -7 K Jk ( -7 hi 2 ) 


( 4 ) 


where Jb, J 7 and Jk denote the spin of the particles B, 7 j and Kj. The notation is close 
to the main application of this paper but we emphasise that at this point the methodology is 
completely general. Assuming the decay to be a series of sequential 1 —>■ 2 decays the amplitude 
can be written in terms of a product of 1 -7 2 HAs times the corresponding Wigner functions 


A(Qb, Qi , ^k\^b, A Ki, Aa' 2 i Ai, A2) ~ 

E Xi<i ( SJ k ) K AkiiJkj D x j ; v (!!,)£ a , iA<2 


( 5 ) 
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where the A j are helicity indices, and 


A^ = Ai — A 2 , (6) 

is a shorthand that we use frequently throughout the paper. The HAs K , K, and C correspond 
to the transitions Bj b —> I\j K r )j , Kj k —> K 1 K 2 and 7 j 7 —> i \£2 respectively. The helicities of 
the internal particles 7 j and Kj have to be coherently summed over. The Wigner D-functions 

D L',m = {jm!\e~ iaJz e~ l ^ Jy e~ i ' yJz \jm) (7) 

are irreducible S'0(3)-representations of dimension 2 j + 1 . The Ji are the generators of an¬ 
gular momentum, and the states |j, m) carry angular momentum j and helicity m and are 
orthonornralised (j, m\f,m!) = 5jj'5 mrn '. To avoid proliferation of indices we denote complex 
conjugation by a bar instead of the more standard asterisk. 

Adaptation to Jb = 0 and K\ = K and K 2 = ir In order to ease the notation slightly 
we move straight to the case B Kj(—> Kir)jj 7 (-^- iih)- 2 The relation D ^~= 
3o,a 7 -a k implies equality of helicities 

A = A 7 = \k ■ (8) 

One may therefore reduce 'H\ 1 \ K —>■ H\, which is the quantity known as the HA in the 
B —> K *£ + £~-literature and carries the non-trivial dynamic information. The HA K.\ k x K 

2 We choose B —> transitions as the main template for the results in this paper. Such transitions are 

governed by the b —» s Hamiltonian, which is the standard in the theory literature and is used to define the 
Wilson coefficients. I 11 the more conceptual sections we shall refer to B —> ATifTtransitions. 
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reduces to a scalar constant (denoted by (JKjKtt) since K\ —> K, K 2 —> n are both scalar 
particles. The third HA £\iA 2 depends on the interaction vertex of the leptons, but is trivial 
to calculate once the interaction is known. We may rewrite the amplitude (5) as 

Jk 

A(B -+ K Jk (^ Kn)^ l x I 2 )) ~ D J x x(n K )D J x XmA J x \ M , (9) 

\=—J K 


where the angles, depicted in figure 1, are Hk = (O,0x 5 O) and 0.£ = ((j)£,6e,— cfre). Note, the 
passage from D to D-functions from (5) to (9) is related to passing from B to B. 

In the lepton-pair factorisation approximation, defined more explicitly in the following sec¬ 
tion, the amplitude A J XXi x ~ -Ha£aiA 2 |j 7 is the product of the hadronic and leptonic matrix 
elements. The angle 4>i is the helicity angle and is usually called simply (j>. Before commenting 
on different conventions of the angles we quote the fourfold differential decay 


d 4 T 


dq 2 d cos 9id cos O^dcj) 


1/2 


ElA 

A1A2 


( 10 ) 


E E A i\,xA'MM D t« («*■) A/d ( q k) Ay m Ay, m . w 

Ai,A 2 =—1/2 A,A'=— J~f 




in terms of amplitudes and Wigner D-functions. For the angles we use the B —> K*££ decay 
as a reference and use the same conventions as the LHCb collaboration [31] (appendix A), 
which differ from those used by the theory community. More precise statements, including a 
conversion diagram, can be found in appendix C.2. 


2.2. Effective Theories rewritten as a coherent Sum of Sequential Decays 

In this section we give the formal steps to derive the expression of the angular distributions. 
The reader interested in the final result can directly proceed to section 3. 

The amplitude (9) is of a completely general form for the decay where yj 7 is an actual 
particle of spin J 7 . In B —^ K*(—t K-k )£\£2 a part of the amplitude is in this form where the 
photon corresponds to the intermediate state (71 = 7 ). In general there are effective vertices, 
so-called contact terms, where the intermediate particles are not present. In the interest of 
clarity we quote the effective Hamiltonian for b —>• sit? 

H eS = c H H eS , t b , H cS = J2 (CiOi + ClO '). (12) 

^ i=V,A,S,P,T 

Above Gf is Fermi’s constant, a the fine structure constant, Ed are Cabibbo-Kobayashi- 
Maskawa (CKM) elements and the operators are 

0 S (p) = s L b £( 75 )£ , Oy(A) = s L ^b £ 7^(75 )£ , O t = s L a^ u b ia^i , (13) 

where O' = O the labels refer to the lepton interaction vertex, qlji = 1/2(1 T 75 ) 9 , 
£,£ -A I 1 A 2 for different lepton flavours and a few additional relevant remarks deferred to ap¬ 
pendix A.2. In passing we add that the notation Og( 10 ) = OytA) is more common throughout 

3 The adaptation from St — > St 2 is trivial and will not be spelled out explicitly. 
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Figure 1: Decay geometries for B —► K*£\l 2 (above) and B —► K*t\£ 2 (below). In both cases i\ = £~, 
t 2 = £~ denote the negatively charged lepton. The conventions are the same as used by the 
LHCb collaboration in [31] (cf. appendix A therein). Comparison to the convention used by 
the theory community can be found in appendix C.2. The pictures are slightly misleading 
in that the angles Qi^k are drawn in the restframe of the lepton-pair and the A'*-meson. For 
decays which are not self-tagging, such as B S ,B S —>• </>(—>• K + K~)[i + ^ at the LHCb, the 
angles (0£, 9k , </>) —t {ir — 9 z,tt — 6k, 2n — <j>), and one can only measure the sum of both decay 
rates (second equation(C.6)). 


the literature. In the case where electroweak corrections are neglected at the matrix element 
level one may factorise the hadronic from the leptonic part. We refer to this as the lepton-pair 
factorisation approximation (LFA) (.B-particle factorisation approximation in the introduc¬ 
tion). Schematically (12) is written as a product of a hadronic part H and a leptonic part £ 
with a certain number of Lorentz contractions between them: 

N 0 Ah N 2 

h cS ~ £ H a c a + ]T + E ■ ( 14 ) 

a— 1 6=1 c=l 

The sum over a, b and c extends over operators with 0, 1 and 2 Lorentz contractions between 
quark and lepton operators. In the example of CyOy = H^O 1 we would have H^ L = CySL^^b 
and = t^l. On a formal level we might think of Oy(Og) as originating from integrating 
out a vector and a scalar particle, in the sense that the Lorentz contraction over index fi can 
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Figure 2: A diagrammatic interpretation of the process, equation (19), used to set up the formalism. 

The decay to two leptons is treated as being mediated by an effective particle r )j i of spin J 7 . 
The factor gKjKi r has no dependence on helicities and depends only on the dynamics of the 
K* decay. 


be written as the sum of products of a spin-one and a timelike spin-0 polarisation vector. This 
is expressed by the well-known completeness relation (e.g. [18,23,30]) 

<T = £ ^(X)Q V (X!)G xy , G\y = diag(l, —1, —1, —1) , (15) 

A,A'e{t,±,o} 

where the first entry in G\x refers to A = A ’ = t and an explicit parametrisation is given by 
w M (±) = (0,±l,i,0)/\/2 , w M (0) = (q z ,0,0,qo)/\/(? , ^(t) = (q 0 ,0,0,q z )/y/(? , (16) 

which is consistent with the parameterisation q^ = (qo,0,0,q z ). The polarisation vectors 
cu(=t,0) are compatible with the Jacob-Wick phase convention [1] (cf. appendix B and the 
corresponding footnote for further remarks). Let us pause a moment and emphasise that 
intermediate results do depend on the convention, which enters the definition of the HAs, 
and this dependence has to be taken into account when comparing to HAs appearing in the 
literature. We choose the convention in [18], since it is compatible with the Condon-Shortly 
convention that is standard for Clebsch-Gordon coefficients and Wigner matrices (e.g. [32]). 

We may think of u as being associated with the Lorentz group 50(3,1). In the rest frame 
q z = 0 the timelike polarisation tensor transforms as a scalar under the restriction of 50(3,1) 
to spatial rotations 50(3). 4 For an effective operator with n Lorentz indices the relation 
(15) can be inserted n times to obtain a HA with n helicity indices. More precisely, the 
direct product of 50(3,1) polarisation tensors decomposes into irreducible representations of 
50(3) polarisation tensors of spin j = 0,..., n and helicities A = —j, . .. ,j. Using the 

expressions in Eqs. A.8 and A.9 the analogue of in (9) on each spin component can 


4 Formally the branching rule for the Lorentz four vector (1/2,1/2) is (1/2, 1/2)so(3,i)|so(3) — 1 (1 + S)so( 3 )- 
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be written as 5 * 

+ J 7 = 0 

+ (HD(Cf)e 1 ^ ■ e^ J 7 = 1 (17) 

(HD{Cf)e 2 ^-e 2 $ J 7 = 2 

where summation over Lorentz indices and the number of operators in (14) are both implied, 
the scalar product is detailed in (A.8) and 

(H^) = (Kj(\)\H^\B) , (£?-»') = (7i(A 1 )£ 2 (A 2 )|£^-^|0) , (18) 


A J ~i — 
•^A,Ai,A 2 “ 


{H a ){C a ) + 


/ ejM\ / o 

\^b)\^'b) e t l e « 


(K)(c 


i A 

-fl to 


are the leptonic and hadronic matrix elements. The helicities in (17) are the helicities of the 
outgoing particles of the HAs, with A for Kj( A) in H b ^ Kj k and \i = \\ — A 2 for ^(Ai)^^) 
in 1 2 . This is the main idea of the formalism: the angular dependence from the ingoing 

to outgoing particle is governed by the Wigner D-function, e.g. e Je ’ x = Z)^ 7 ^ (fi#) for 

£ 7 j 7 (A)-»H(Ai)t 2 (A 2 )^ w j 1 j c j 1 j s inherent in (3). The generalised HA then becomes essentially 
a sum over all spin components J 7 necessary to saturate the Lorentz indices in the effective 
Hamiltonian, 

/ 2 j n min(J 7 ,J K ) 

-> ic, K {-+K*)iJ 2 ) = v * + ^ ^ 19) 

J 7 =0 A=—min( J 7 , J/^) 


where the overall factor follows from (3). A schematic representation of equation (19) is given 
in figure 2. The differential decay distribution (10) is replaced by a similar expression 


d 4 T 


dq 2 d cos Oid cos Oxdcj) 


DA 


AiA 2 

( n *> («*■) D i\ (^> D i,x, m . 


2 Jk + 1 
47T 


EEE 

A 1 A 2 J'y A J^A' 


with additional coherent sums over the spins J, 
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1/2 2 min(J 7 ,JKj) 

E- E ■ E=E E 

AiA 2 Ai,A 2 =—1/2 J 7 A J 7 =0 A=— min( J 7 , ) 


( 20 ) 


( 21 ) 


and likewise for the sum over J’. Ah 


3. Angular distribution and Wigner .D-functions 

We now apply the method introduced in the previous section to decays governed by the b —> 
si \(-2 effective Hamiltonian (12). First we consider the decay B —> K* (—> Kir) l\t 2 , and then in 
section 3.3 we present similar results for B —> Kl\l2 ■ The related decay Aj, —> A (— > Nir) i\£.2i 
where N = ( p,n ), can also be treated within this formalism, and will be briefly considered in 
appendix E. 

5 In the notation used throughout the literature H l = (H^e 0 ^ 0 = is known as the timelike HA [23,30]. 

By virtue of the equation of motion the timelike HAs can be absorbed into the scalar and pseudoscalar HAs, 

cf. appendix C.5. 
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T S[P] 

py[4] 

y t ip] 

r x 

MX) 

14 ) 75 ] 

0 

1 

,1 Ax r, iAxi 
<7 UJiu/ l 

1 


Table 1: The definitions of the T A and their associated spin J 1 {X). The contributions J 7 (JT) = 0,1 
give rise to the Si- and Pi-wave amplitudes respectively. The basic polarisation vector 
is given in (16) and the composed ones can be found in equation (A.5). The precise value 
of the helicity index Ay is specified when the leptonic and hadronic HAs are defined in 
Eqs. (A.13,C.14,D.4). Note that the additional structure r Ts = 07 , 1,75 can be absorbed into 
the other tensor structures due to the identity 0" a/3 75 = — | e al3 ^ t '(with the £0123 = +1 
convention for the Levi-Civita tensor). Timelike contributions 7 ^ [75 ]w M (f) can be absorbed 
into T s ’ p respectively, as detailed in appendix C.5. Above = */ 2 [ 7 ^, 7 ^]. 


3.1. B -A K* (—>■ Kir) 44 

The use of the effective Hamiltonian (12) in the LFA restricts the partial waves to J 7 = 0,1 
terms in equation (17). The discussion of higher partial waves ( J 7 > 2) is deferred to section 5. 
The matrix element for (12) is then given by the sum of an Si- and Pi -wave amplitude (with 
the subscript £ referring to the partial wave in the angle 81 ): 

Mx u a 2 = (K*(-a Kir)h{\i)U^2)\H eS \B) 


V3 

47T 


,a 2 4),o (44 8. 


U X 4- 




( 22 ) 


where the hat denotes the effective Hamiltonian without the ch prefactor (12). There is no 
D -wave since the two-indices in the tensor operator (12) are antisymmetric and therefore in 
a spin 1 representation (cf. discussion in section 5.2 on higher spin operators). The K* has 
spin 1 and is therefore always in a Pk -wave in the 0jy-angle, with analogous meaning for the 
K subscript as before. Above we have used Hq X( (12) = Sox e to impose 4 iA 2 on the scalar part 

of the matrix element. The principal objects to be calculated are the amplitudes A XXl x 2 - For 
H eS (12) the Si- and Pi-wave amplitudes (that is to say A 0 and A 1 respectively) are written 
as 


A 0 

A),Ai,A 2 

A 1 

AMiAz 


ttS rS 

V r v 


+ H P C P 


-me 


A '°Ai,A 2 


-me 


Ai,A 2 

A r A 


A A -'Ai ,A 2 


+ H T x e T XlM 


2 H^e 


T t 

Ai,A2 ’ 


(23) 


with the relative signs and factor of 2 emerging from the (double) completeness relation (A.3), 
and the leptonic and the hadronic HAs are 

H? = (K*(X)\sT X b\B) , e Xl \ 2 = (h(\i%(^W x £\0) , (24) 


the expressions in (18) contracted with the corresponding polarisation vectors. Explicit results, 
as well as a more precise prescription concerning r4 are given in Appendices A.3 and C.5 
in Eqs. (A. 13) and (C.14) respectively. Squaring the matrix element in (22), summing over 
external helicities and averaging over final-state spins, one obtains an angular distribution 


Ik* ( q 2 , Mr, 4) 


32vr d 4 r 

3 dq 2 dcosOi dcos0K d<f> 


32t r 

ir 


■N X] |A^Ai,A 2 

Ai,A 2 


2 


(25) 
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with Ik* being a shorthand and 327 t/ 3 is a convenient normalisation factor. The factor Al, 


AT = 


CH | ^kin 


^kin — 


2 6 7^ 3 m^q ,2 


(26) 


is the product of the prefactor resulting from the effective Hamiltonian ch (12) and the kine¬ 
matic phase space factor. The matrix element is defined in (22). Above A b = \(m 2 B , m 2 K *, q 2 ) 
and A 7 * = \{q 2 , m 2 ^, ) where A (a,b,c) is the Kallen-function defined in (B.l) and related 

to the absolute value of the three-momentum of the K* and the lepton pair by (B.2). 


3.1.1. Angular Distribution 

The squared matrix element initially contains a plethora of different products of four Wigner 
functions. However, these correspond to pairs of direct products that can be reduced to single 
Wigner functions by the Clebsch-Gordan series 

j+i J J 

J2 E E C Mn,p C Nm D il,N^) ■ (27) 

J=\j-l\ M=—J N=—J 


Applied separately over the angles II k = (0, 6 k, 0) and Up = (</>, &e, —(f>), along with the identity 
D l m m' (4 = (—l) m ~ m D l _ m _ rn , (H), this allows the angular distribution to be written in the 
compact form 


4 °* {q 2 ^kMc) 


Re 


n 0,0/ 2'io0,0 

Urn 


V)4 


Gl’°(q 2 n 


Gf{q 2 n 


2,0 

0 

2,2 


T^oV^o’ 1 

+ Gf(q 2 )I % 1 

+ 4’V)4’ 2 


+ G° 0 ’ 2 (q 2 )Q °’ 2 

+ 4’V)4’ 1 

+ G 2 /(q 2 )n 2 / 


+ 

+ 


(28) 


where the superscript (0) is a reminder that only Sg- and Pa -wave contributions were used to 
describe the amplitude (22). The angular functions H are given in terms of Wigner D functions 

44 = 44 (n K ,n e ) = d 1 « >0 (n K ) < 0 (n e ) = d% i0 (4) 4U ( 4 ) . ( 29 ) 

The variables Q' K = (4>, 6k, —0) and 4 = (0, 6e, 0) form an angular reparametrisation that will 
prove convenient when we discuss partial moments. The label Ik corresponds to the (Ktt)- 
system, lg to the dilepton system, and the common index m is the azimuthal component <f> of 
either partial wave. The observables Gm’ le are functions of q 2 and the relation to the standard 
observables in the literature is given in section 3.2. The explicit Wigner D -functions used 
above are given by 

4, 0 (fi) = 1 , 4,0 (4 = \(1 cos 2 6-1) , 4,o (4 = vf e ' 2 ' 0 sin2 6 ’ 

4o(4 = cos0, 4o(4 =—^=e“^sin0, 4o(4 =-\/|e _< *sin20 , (30) 

V 2 ’ V 8 

and can be related to spherical harmonics Y[ m (6 , (f>) or associated Legendre polynomials Pi m (x) 
as 

4n,o 4 0, -4>) = \f^-Y lm (6,4>) = ^^PUcos 6)e~ im * . (31) 


9 









We comment briefly on four features of the angular distribution (28), all of which are encoded 
by the double Clebsch-Gordan series (27) but which can also be seen to emerge from the 
underlying physics: 

• The second helicity index of all Wigner D-functions in the angular distribution is zero. 
The latter is the difference of the helicities of the final-state particles, which is zero since 
these helicities are summed incoherently, (Ai — A 2 ) — (Ai — A 2 ) = 0. 

• The first helicity index m is identical in all pairs of Wigner D-functions appearing in the 
angular distribution. This index contains the helicities of the internal particles, summed 
coherently. One can also see this as a property of the freedom of defining the reference 
plane for the angle cj>. 

• The range of the indices Ik and l( is fixed between the range 0,..., 2 max[J^]. Including 
only J 7 < 1 contributions emerging from the dimension-six effective Hamiltonian (12) 
hence imposes 0 < lg < 2, and likewise Jr = 1 imposing 0 < Ik < 2. 

• The absence of angular structures with Ik = 1 is specific to this decay, due to the final 
state consisting of (pseudo)scalar mesons. 

The first three features are universal to such decay chains and apply even if some of the particles 
involved are fermions, for example in the decay A& —> A (-> (p. n)n) £±£ 2 , see appendix E. 


3.2. Relation of the G lK,tl to standard literature observables 

rri 

The functions Gm’ li , omitting the explicit ^-dependence hereafter, are defined in terms of the 
standard basis of observables 9 i{q 2 ) parametrised in (C.l) by 


Co ° — g (3 (gic + 2gis) — ( g 2 c + 252s)) > Gq' 1 — - (g§ c + 2g§ s ) , 
Co 0 = ^ (6 (jjic - gis) ~ 2 ( g 2 c - g2s)) , Go" 1 = ^ (gee - ges) , 

r 2,i _ }6_ r r 2,2 _ 32 r 

— yghs j — 3 h4 , 


Gq’ 2 — — (<72c + 2<72s) j 

99 32 

Gf=-{ 9 2c- 9 2s) , 

99 32 

Gf = y Gs , (32) 


where we have defined £ 3 , 4,5 = ( 53 , 4,5 + * 59 , 8 , 7 )- 

The twelve quantities (32), keeping in mind that the last three are complex, have been 
rewritten in several ways in the literature. A frequently-used form is the set of observables 
given in [33], constructed to be insensitive to form factors. In the notation of LHCb [26], which 
includes their, and therefore our, angular conventions, the observables are given in terms of 
Gm' i( by: 6 


(-Pl)bin 

( P 3)bin 


^Re 

G 2 / 

) bin 

LHCb 

■A/bin 


l Im 

G 2 / 

) 

\ 


/ bin 


LHCb 2A/bi n 


(P 2 ) 


bin 


LHCb 


2Gj’ 1 - Gj 1 } 

_ / bin 

3A/bi n 


6 The extension of these relations to CP-odd and CP-even combinations, in the spirit of [34], is straightforward, 
see section 4 of [33]. 
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(Pi) 


bin 


Re 


Gt 


2,2 


bin 


<^)bin 

where we defined 


LHCb 


LHCb 


Re 


^bin 

gI’ 1 


2V2.Nl 


u 


bin / p'\ 

’ V 6/bin 


Im 


Gt 


2,2 


bin 


LHCb 


LHCb 


^bin 


Im 


Gt 


2,1 


bin 


2V2.N1 


bin 


( 33 ) 


W)) hin = dVf(V), 

J bin 

as the integral over q 2 bins of the observable of interest, and' 

Nbm = 4 ( G°v - ; G ; 


1 


VV • -^bin = \l~( a Y - 


bin 


2,2 


> 


bin 


G^ + G, 


2,2^ 


bin 


(34) 


Three other combinations of the Gm ,le can be related to the branching fraction jjV , the forward- 
backward asymmetry A^b and the longitudinal polarisation fraction Fj j [35]: 


dT 

dq 2 


bin 


( F l) bin = 


_ 3 / „ 0 ,o\ 

“ 4 \ G ° /bin 




(^fb) 


bin 


1 \°iib te 

or) 


0.1 \ 


LHCb 2 /,o0,0\ 


bin 


3 (G°’° 


(35) 


bin 


The observables in Eqs. (33,34,35) correspond to the twelve gi. The definitions of the P[ above 
correspond to those used by LHCb [26]; we give the correspondence to the observables defined 
in [33] in appendix C.2. 


3.3. B KIV 2 

Having shown the B —> K*£ 1^2 HA analysis in detail we are going to be rather brief on 
B —> K£ 1^2 ■ Skipping the step in (5) we directly write down the Si- and Pi -wave amplitudes 
(analogue of equation (23)): 


A 0 

"^0,Ai,A 2 


S r s 


h a C 


Ai,A 2 


+ h p C p 


h v r v 


Ai,A 2 > 

h A r A 
~ n 


+ ^ T ^Ai,A 2 


— 2 h Tt P Tt 
m ^'A 1 ,a 2 


(36) 


where the C* Ar) are the same as in the B —> K*£\l 2 decay, and the hadronic HAs are taken 
over the same set of operators, but defined instead for B —>• K transitions. We again refer the 
reader to appendix A.l for a clarification of the signs and factor of 2 that emerge from the 
(double) completeness relation. 

The reduced matrix element is then the sum of the Si- and Pi -wave amplitude 


M 


Ai,A 2 


V~Vk 


(”^0.Ai,A 2 ^A 1 A 2 + "4o,Ai,A 2 Ai,A« (^)) 


7 In terms of the g, (g 2 ) basis, W bi „ = ^ (ff 2 s) bin and AC, = fV~ ( 92 c) bin (ff 2 S ) bin . 


(37) 


11 



























where fig = (0, 0g, 0) in this case. The angular distribution (with 0 < 9g < it) is given by 
squaring the matrix element 


Ixiq 2 ,9g) 


d 2 T 

dq 2 dcosOg 


* E 

Ai,A2 


(38) 


Using (37) one obtains 

4°) = g(°4 2 ) + m + G^(q 2 )D 2 0fi (fig) 

= G(°\q 2 ) + G ( ' 1 \q 2 )Pi(cos9g) + G^ 2 \q 2 )P 2 ( cos9g) 

= G^(q 2 ) + G^(q 2 ) cos 9g + G^(q 2 )^ (3 cos 2 9 t - l) , (39) 


where we used Pi( cos Qg) = D l 0 0 (fig) and ZDq 0 (fig) = 1. For convenience, we have given results 
in terms of the explicit angle 9g using equation (30). The superscript (0) is again a reminder 
that the restriction to lg < 2 is a consequence of only including Sg- and Pg- waves in (37). The 
explicit functions G^ 0,1,2 \ whose ^-dependence we omit hereafter, are given in appendix D in 
equation (D.2) in terms of HAs. 8 

With respect to the parametrisation of the angular distribution used in the experimental 
community, [37] 


i dr 

r d cos 9g 


- (1 - F h ) (1 - cos 2 i 


+ 2^ h ^fb cos( 


the relationship to the G® in (39) is given by 


T = 2<G(°)> , A fb 


(g (1) ) 

• Y 2(G(°)) ’ 


Fh 


(G(°>) + (G&) 
(G(°)) 


(40) 


(41) 


where ( X ) = J dq 2 X denotes the integration or appropriate binning over q 2 and a = ±1 
depending on the conventions. 9 


4. Method of total and partial Moments 

The MoM is a powerful tool to extract the angular observables Gm’ le by the use of orthogonality 
relations. In B physics, for example, the method has been applied to B —> J/T (—> ££)K*(—> 
Ktt) type decays [29] during the first B-factory era. 

In experiment the angular information on B —> K*££ has been extracted through the like¬ 
lihood fit method, at the level of 4* [27], and it has also been suggested for analysis at the 
amplitude level [38]. A possible advantage of the MoM over the likelihood fit is that it is 


8 The observables G (i d and the angular coefficients used in the literature [36] are related by a(q 2 ) = G^ — \G^ 2 \ 
b(q 2 ) = G^ and c(q 2 ) = §G (2 ^ where = a 4- bcosdi + ccos 2 6e. 

9 In our conventions by definition ctghz = 1 and the translation to the LHCb conventions [37] are as follows 
(tghz = a(B ± ) and ctghz = —a(B°,B°). The charged and neutral decays are different because the neutral 
mode, being observed in Ks, is not self-tagging. Comparing with the theory paper [36] we find ctghz = 
—obhp f° r both charged and neutral modes. 
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less sensitive to theoretical assumptions. More precisely, one can test each angular term inde¬ 
pendent of the rest of the distribution. Generically the fourfold angular distribution can be 
expanded over the complete set of functions f im’ li (29) 


Ik* (q 2 ,QK,^e) = ^ 


min (IkM) 

E Re 


0 m =0 




(42) 


of which the distribution I^l (28) is a subset. Note that the sum over m does not need to be 
continued for negative values since Ik* is real-valued. By using the orthogonality properties 
of the Wigner D-functions (e.g. [39]) with fl = (a,/3, 7 ) 

/ I /*27r /*27 t o 2 

^dcosfi Jv da d 7 D 3 m n (Q) D l p q (fi) = + 1 6 j i6 mp 6 nq , (43) 

the MoM allows to extract the observables Gm’ le from the angular distribution. In particular 
one can test for the absence of all higher moments and therefore test very specifically the 
assumptions made when deriving the distribution I^l (28). We refer to this method as the 
method of (total) moments or simply MoM with results given in section 4.1. Integrating over 
a subset of angles, referred to as partial moments, is discussed in section 4.2. In the latter case 
orthogonality does not hold in the generic case and different Gm’ le enter the same moment. 

Elements of the MoM have previously been applied to Aj, —> A (—» ( p , n) 7 r) t\l 2 [40] and 
more systematically to the other channels discussed in this paper, crucially including a study 
of how to account for detector-resolution acceptance effects, in [24]. Our study differs from 
the latter in that we start at the level of the HAs, and obtain the distribution (42) through 
a direct computation, whereas the other studies proceed backwards and directly expand the 
decay distribution in the orthogonal basis of associated Legendre polynomials. Our approach 
is therefore advantageous in that it provides additional insight, by clarifying the structure of 
the decay distribution (28) and what type of physics goes beyond it. This is an aspect we 
return to in section 5. 


4.1. Method of total Moments 

In order to condense the notation slightly we define the scalar product 

1 /*1 pi p 2 tt 

(f(SI)\g{ty)e K Ot<l> = ^ J d cos 6 k J d cos J d<j)f(n)g(n), (44) 

normalised such that (1|1) = 1. Using (f(Q)\g(Q))g K o ec f ) we can thus extract all observables 
Gm’ 1 separately from each other, by taking moments 10 

M l *b = (n l «’ l ‘\I K *(q 2 ,n K ,n e ))o K e e t = C 1 ^G 1 ^ , (45) 

where 


Jx,h — 


1 + dmO 


2{2l K + l)(2k + l) ' (46) 

Using the equation above the terms in (28) are given in Tab. 2. Furthermore, the orthogonality 

10 The moments M l m’ e and the quantities Si e ,i K: m introduced in [24] are related as follows: STrGg’ 0 ^^^^, m = 

k _ T\/r l K^e. 

lvl m / I'm 
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M 1k ’ lt 
±v±m 

a 

oo 

o 

o' o 

o 

rH|CO 

1 /^ 0,2 
5^0 

1 /^# 2,0 
5 b 0 

i n 2 ’ 1 
15 

1 ^ 2,2 
25^0 

i r 2 ’ 1 
30^1 

1 /^ i 2,2 

50^1 

1 K ' 2,2 
50^2 


Table 2: Moments in terms of G\%' le as defined by equation (45) with factor of proportionality 

c[%’ le evaluated with (46). 


condition also implies that 

M™' = 0 , Vm and j > 3 or j' > 3 , 

M l J= 0, V/,m. (47) 

Hence the higher and Ik = 1 moments vanish, providing a very specific test of the theoretical 
assumptions behind 1 ^ 1 . 


4.2. Partial Moments 

The results given previously show how to extract the individual Gm’ le ■ We propose the method 
of partial moments whereby one integrates only over a subset of angles. The distributions might 
be regarded as generalisations of uni- and double-angular distributions as these in turn can be 
viewed as partial moments with respect to unity. The method is effectively a hybrid between 
the likelihood fit and the total MoM. To this end we define the further scalar products 

(f(Q)\g(n)) 0 ( l> = — J ^dcos 0 J ^ d(j)f{£l)g(n), 

(f(n)\g(n)) eK 0 e = | j dcos9 K J dcos 0 ef(Sl)g(n) , (48) 

again normalised such that (1|1) = 1. The orthogonality relation (43) can then be rewritten 
as 

{D pfi (^) \^m,0 (^())0e4> = 21 _|_ ^jrfrnp ■ (49) 




4.2.1. Integrating over 6^, 4>: 

The partial moment over 9i and </> is defined and given by 

*4(M = (Di, o m I Ik- n K , «<) W = ^Tl) S D ™fi (^) 

y i K >o 

Assuming the distribution (28) (Ik = 0, 2) there are six non-vanishing moments 
k° 0 (e K ) = G°’° + G 2 ’°L> 2 o (n K ) = G°’° + \ (3 cos 2 9 K - l) G 2 ’ 0 , 
kl(9 K ) = \ (C 0 ,’ 1 + G 2 /Dl 0 (Q k )) = J + \ (.3 cos 2 0 K - l) G 2 ’ 1 ) 

kl(9 K ) = ^ (g°’ 2 + GfDl o (fi^)) = J (g°’ 2 + \ (.3 cos 2 9 K - l) G 2 ’ 2 ) 

1 ' 

-i 
6 


k\(9 K ) = \G 2 { X D\^K) 


— 1 /3 . „ „ 2,1 

— \ - sm 29kG i , 

6 V 8 11 


(50) 
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( 51 ) 


k 2 i(e K ) = ^Gl’ 2 Dl fi (n K ) 
kl{e K ) = ^GfDl fi (Sl K ) 


~ ■*" r < 2 > 2 

10 V 8 1 

1 /3 . 2 .-,2,2 

= —i -sin Ok G, . 

10 V 8 2 ’ 


where we used Dg 0 (^at) = 1- As was the case in the MoM, with respect to the distribution 
I^l higher partial moments vanish 

tfn(@ K ) = 0 VZ*>3 ,Vm. (52) 


4.2.2. Integrating over 9k, <$>'• ^(^)-nioments 

The partial moment over Ok and <j) is defined in complete analogy with the previous partial 
moment (50) by, 


tem = {d% o w i iK-tfAMVox* = 


1 + 5m0 

2 (2 l K + 1) 


E D m,o (Uk)#*’ 1 *, (53) 


h>o 


where we make use of the reparametrisation of angles given in (29). Again assuming the 
distribution (28) (l £ = 0,1,2) there are four non-vanishing moments 

l° o (0 £ ) = G°’° + G^Dl o (fij) + G° 0 ’ 2 Dl o (nj) 

= G°’° + cos 0,Gg 1 + i (3 cos 2 0, - 1) G°’ 2 , 

*o(fc) = ^ (g 2 0 ’° + Gq 1 Ho 0 (ni) + G 2 q 2 Dq Q (n£)) 

= * (g 2 ’° + cos 0 £ G 2 + + i (3 cos 2 0 £ - 1) G 2 ’ 2 ) , 

ilm = ^ (g 2 ’Xo (nj) + G 2 { 2 Dl o (O')) 


-i 


sin 0,G 2,1 + \/ ^ sin 20,G^ 


2,2 


10^2 

<?(»<) = ^G} 2 Dl„ (SJJ) = 4y/| S in 2 o 


2,2 


(G 2 , 


(54) 


where we used Hg 0 (O^) = 1. With respect to the distribution higher partial moments 
vanish 

= 0 , V/jf > 3, Vm and Ik = 1, Vm . (55) 


4.2.3. Integrating over 6k, Of- m °ments 

Finally, we can consider projecting on to moments of the form D l m 0 (fix) D 1 ^, 0 (O^) with 
respect to Ok, 0 £ . In this case the full orthogonality relation (43) no longer holds, but due to 
(27) there exist selection rules as to which of the Gm’ le can contribute to the partial moments 

p'XlZ'W = (D l «,o (0,0 k,0)D 1 ^, o (0,0,, 0) I lK*tf,n K ,n t )) eK e t . (56) 
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Assuming / 


(0) 


K* 


a few non-vanishing moments are 


Po’fiW) = (<’° + Refe - ^^’ 2 ]) , p{j;£(0) 

Poflit) = ^ (< 2 - Re[e~ 2l4 G 2 /}^) , p 2 o ' O o (0) 

Poi^) = ^, Po 2 0 (^) 

Pifa) = e[e~ lct> G 2 /] , pl’fa) 


1 r , 0,1 
3 G ° ’ 

^ (6 G 2 0 ’° - R e[e~ 2i +Gl' 2 ]) , 

+ Re[e- 2 ^’ 2 ]) , 
iRe[e-^G?’ 2 ] . (57) 


A consequence of the fact that the full orthogonality of the Wigner functions has been lost is 
that higher moments contain lower G- functions. As an interesting example we quote 



1 

9\/l0 




9v/10 56c ' 


(58) 


This quantity is of some interest since g§ c = 0 in the SM, as it involves scalar and tensor 
operators at the level of the dimension-six effective Hamiltonian (12). 


5. Including higher Partial Waves 

The compact form of the angular distribution I^l (28) is a consequence of the LFA and 
the restriction to the Pk -wave in the (it7r)-channel. In this section we elaborate on the 
consequences of going beyond these approximations. The double partial wave expansion is 
outlined in section 5.1 followed by a qualitative discussion of the effect of higher spin operators 
and the inclusion of electroweak effects in sections 5.2 and 5.3 respectively. In section 5.4 we 
emphasise how testing for higher moments can be used to diagnose the size of QED corrections. 
Throughout this section we change the notation from £\£ 2 —> £ + £~ for the sake of familiarity 
and simplicity. 

5.1. Double Partial Wave Expansion 

In order to discuss the origin of generic terms in the full distribution (42), it is advantageous 
to return to the amplitude level. Somewhat symbolically we may rewrite the amplitude (19), 
omitting the sum over J 7 , as 

A(B Kj(A)(-> Ktt)£+(A 1 )£-(A 2 )) = A{^ K D j « q (n K ) D) ^ (Q,) (59) 

with Ai = A\ — A 2 as defined in (6). The two opening angles 9k and allow for two separate 
partial wave expansions. The partial waves in the 9k- and ^-angles are denoted by Sk, Pk , • ■ • 
and Si, Pi,... respectively. 

Throughout this work we mostly restricted ourselves to Kj = K* thereby imposing Jk = 1 
i.e. a Pk- wave. The signal of K* is part of the ( Kk ) Pk- wave. The importance of considering 
the S^-wave interference through J\g(800) (also known as k( 800)) was emphasised a few years 
ago in [41]. The separation of the various partial waves in the (A"7r)-channel is a problem that 
can be solved experimentally e.g. [42], We refer the reader to Ref. [19] for a generic study of 
the lowest partial waves at high q 2 . 
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Figure 3: Examples of virtual QED corrections to B —► K£ + £~, where either a photon is exchanged 
between the decaying &-quark and a final state lepton, with effective operators Oy,A (left) 
and 0 7 (middle), or a second photon is emitted by the charm loop (right). Other topologies 
relevant for higher moments include the interaction of the leptons with the spectator as well 
as the B- and A’W-meson. 


The second partial wave expansion originates from the lepton angle 9g, which will be the 
main focus hereafter. By restricting ourselves to the dimension-six effective Hamiltonian equa¬ 
tion (12) as well as the lepton-pair facorisation approximation (LFA ) 11 only Sg- and P^-waves 
were allowed (cf. equation (22)), bounding lg < 2 in (42). This pattern is broken by the inclu¬ 
sion of higher spin operators and non-factorisable corrections between the lepton pair and the 
quarks. It is therefore important to be able to distinguish these two effects from each other. 

5.2. Qualitative discussion of Effects of higher Spin Operators in iT eff 

Operators of higher dimension are suppressed and neglected in the standard analysis. Opera¬ 
tors of higher spin in the lepton and quark parts are necessarily of higher dimension and bring 
in new features. An operator of (lepton- and quark-pair) spin j is given by 

c>0) = s L rg-^, b It^ (eo) 

with ' ‘ ‘ = D±D, with D the directional covariant derivative and 

curly brackets denoting symmetrisation in the Lorentz indices. In passing let us note that in this 
notation = Oy = O 9 with Oy defined in ( 12 ). The operator (60) is of dimension d 0 (j) = 
4 + 2j and the corresponding Wilson coefficients are suppressed by powers of mw- Neglecting 
electroweak corrections and including the dimensional estimate of the matrix elements the 
leading relative contributions are given by (rrib/mw) 2 ^ ^ where (rrib/mw) ~ 6 • 10 -3 . 

Operators of the form (60) present new opportunities to test physics beyond the SM provided 
that their contribution is larger than that of the breaking of lepton factorisation through 
electroweak corrections. The operator O <y2 \ for example, gives rise to non-vanishing moments 
of the type G 2’ 1 in B — > K*(—f Ktt)I + £~ and G 4’ 4 in B —> /Q(— Ktt)I + £~ , [44] both of which 
are absent in the LFA. 
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5.3. Qualitative discussion of QED Corrections 


The B —> K£ + £~ channel allows the discussion of the consequences of going beyond the LFA in 
a simplified setup, and is of particular relevance because of a recent LHCb measurement [25]. 

In the LFA (38) the single opening angle 0f of the decay is restricted to If < 2 moments 
in I (42). More precisely, If < 2 j with as in (60) (see also the discussion following 
equation (28)). From the viewpoint of a generic 1 —> 3 decay there is no reason for this 
restriction, as it is only the sum of the total (orbital and spin) angular momentum that is 
conserved. However, in the LFA the B —> I\[£ + £~] decay mimics a 1 —> 2 process, imposing 
this constraint. This pattern is broken by exchanges of photons and W- and Z-bosons, as 
depicted in figure 3 for a few operators relevant to the decay. The W and Z are too heavy to 
impact on the matrix elements, but their effect is included in the Wilson coefficient. 

As stated above QED corrections turn the decay into a true 1 —> 3 process, and this neces¬ 
sitates a reassessment of the kinematics. By crossing the process can be written as a 2 —> 2 
process 

B(p B ) + £-(-£ 1 )^K(p) + £-{£ 2 ) , (61) 

with Mandelstam variables s = (p + Q) 2 , t = (£\ + f^) 2 = T and u = (p + l\) 2 , 


2 


{m 2 B + m 2 K + 2 mj - q 2 ) ± Pi\J\ ( m 2 B ,m 2 K , q 2 ) cos Of 


(62) 


obeying the Mandelstam constraint s + t + u = m 2 B + m 2 K + 2m 2 . Crucially, the kinematic 
variables s and u become explicit functions of the angle Of. In a generic computation these 
variables enter (poly)logarithms, which when expanded give contributions to any order If in 
the Legendre polynomials. This statement applies at the amplitude level and therefore also to 
the decay distribution (39) 


d 2 T(B -> K£+£r ) 
dq 2 dcosOf 

The B —> K££ moments are simply given by 


Y.G^Pi^osOf) . 
h>0 


M lf )= / d cos 6fPig (cos Of) 


d 2 T{B ->• K£+£~) _ 1 

dq 2 dcosOf 21f + 1 


/ -Ah) 

(jr a 


(63) 


(64) 


where we have introduced a lepton-subscript for further reference. In the SM the effects are 
dependent on the lepton mass, for example through logarithms of the form In [mf/m^] times the 
fine structure constant. There are new qualitative features of which we would like to highlight 
the following two: 


• Both vector and axial couplings Oy^ = 09 ( 10 ) (12) contribute to any moment If > 0. In 
the LFA l^-odd terms (forward-backward asymmetric) arise from broken parity through 
interference of Oy and Oa (12). The physical interpretation is that there is a preferred 
direction for charged leptons in the presence of the charged quarks of the decay. In the 
specific diagram figure 3 (left) it is the charge of the 6-quark which attracts or repels the 
charged lepton(s) with definite preference. It is possible that one can establish a higher 
degree of symmetry by using charge-averaged forward-backward asymmetries. 

lx We remind the reader that in the LFA no electroweak gauge bosons are exchanged between the lepton pair 
and other particles when calculating the matrix element. This is the same approximation that is relevant 
for the endpoint relations [18,43]. 
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• A key question is how the moments vary in It- In the absence of a computation a precise 
answer is not possible. Nevertheless we can assess the question semi-quantitatively by 
considering for example the triangle graph between the photon, a lepton and the 6-quark 
in figure 3 (left) and the corresponding one with the s-quark. Neglecting the Dirac 
structures the triangle graph is given by Co(mj,p'^,s[cos9t],mj,0, mf). 1213 Expanding 
this function in partial waves Cq = Xm>o ( cos we find that |Cq^| does fall 

off in It- Averaging over several configurations (cf. footnote 13) we conclude that the 
It = 2 (D- wave) contribution is suppressed by approximately a factor of 2 with respect 
to It = 0, with a slightly steeper fall-off with increasing It for the 6-quark versus s- 
quark vertex correction. Note the graph where the photon couples to the other lepton 
comes with a different Dirac structure and is not obtainable through a straightforward 
symmetry prescription. We therefore think that it is sensible to consider those graphs 
separately. We stress that this semi-quantitative analysis does not replace a complete 
QED computation, which would include corrections to Wilson coefficients, all virtual 
corrections and importantly also the real photon emission. 


We now turn to the most important consideration, the relative size of the QED corrections 
versus higher spin operators. For effective field theories of the type (H eS ) C&{n F ){OW(ji F )) 
(12), the precise separation scale fj ,f is arbitrary to a certain degree and effects are therefore 
contained in the Wilson coefficients as well as the matrix elements. We find it convenient to 
discuss the effect at the level of the Wilson coefficients. For the latter QED corrections arising 
from modes from m\y to ^f — nr b can be absorbed into a tower of the higher spin operators 
(60). The leading contribution to the corresponding Wilson coefficients from the initial 
matching procedure and the mixing due to QED behaves parametrically as 


C® 


0 ( 1 ) , lf 0 ( 1 ) 

‘ \ml) K„,)l 


for j > 1 , 


(65) 


where we have implicitly used /ip = in (H eS ) ■ Above a is the 

fine structure constant and fj parametrises the comparatively moderate fall-off of the higher 
moments due to QED. In the SM one therefore expects QED effects to dominate over those 
due to higher spin operators, except for j = 2 where they could be comparable [44]. At the 
level of matrix elements this hierarchy could even shift further towards QED as a result of 
infrared enhancements through h^m^/mQ-contributions. 

The discussion of B — > K*(—> Ktt)£ + £~ is similar, but involves the kinematics of a 1 —> 4 
decay. The decay distribution becomes a generic function of all three angles 9t, 9 k and 4>. It 
should be added that the selection of the K* —> Ktt signal (P^-wave) restricts Ik = 0,2. 


5.4. On the Importance of testing for higher Moments for B ->■ K^£+£~ 

We have stressed throughout the text that it is of importance to probe for moments that are 

vanishing in the decay distributions I^l (28) of B — > /!*(—>• Kn)l£) and 1$ (38) of B —>• Kll 

12 We use conventions for the Passarino-Veltman function Co(Pi,P 2 iP 3 >mf,ml,such that the two-particle 
cuts begin at pi > (mi + m 2 ) 2 , p\ > (m 2 + m 3) 2 and > (m 3 + mi) 2 . 

12 We have refined this analysis by taking into account that the b- and s-quark only carry a fraction of the 
momentum of the corresponding mesons. This amounts to the substitution p% —> (pg — xp ) 2 and s[u] —> 
(£2 [-h] + xp ) 2 with x being the momentum fraction carried by the s-quark. For the vertex diagrams one 
expects the Feynman mechanism (i.e. x ~ 0) to dominate. This changes when spectator corrections are 
taken into account. 
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respectively. In this section we highlight specific cases of current experimental anomalies in 
exclusive decay modes where their nature might be clarified using an analysis of (higher) 
moments. 


5.4.1. Diagnosing QED background to Rk 


In the SM the decays B + —> K + e + e and B + —> K + p, + p, are identical up to phase-space 
lepton mass effects and electroweak corrections. The observable 

^ ^ _ B(B+^K+p+p~) 

K Itemin’tfmax] B(B + —> I\+e + e~) 


q 2 . , 
’min 


9max] 


( 66 ) 


has been put forward in Ref. [45] as an interesting test of lepton flavour universality (LFU). 
Above i n / max stands for the bin boundaries. Neglecting electroweak corrections the SM pre¬ 
diction is Rk | 6 j Q e y 2 — 1-0003(1) [46], which is at 2.6cr-tension with the LHCb measurement 
at 3fb _1 [25] 

R k = 0.745lg;^(stat) ± 0.036(syst) . (67) 

Previous measurements [47,48], with much larger uncertainties, were found to be consistent 
with the SM as well as (67). This led to investigations of physics beyond the SM with Cg e ^ Cg ^ 
(where O^f = b r y a sl'y a £) amongst other variants for which we quote a few recent works [49-58] 
as well as the general review [59] for further references. 

Let us summarise the aspects of QED corrections which are of relevance for the discussion 
below: i) they break lepton factorisation and therefore give rise to higher moments, and ii) 
they depend on the lepton mass, for example through logarithmic terms of In {mi/mb)- In view 
of the lack of a full QED computation 11 we suggest diagnosing the size of QED corrections, as 
well as their lepton dependence, by experimentally assessing higher moments. 15 The latter is 
directly relevant for Rk- Let us be slightly more concrete and define the normalised angular 
functions as follows = G^/(2Ggj}) (63) (in this convention 2 G^ = dT(B —> Kl€)/dq 2 , 

G^ = Afb and = ( Fh — l)/2 in the notation of [36]). We would like to stress the 
following points: 


• How to distinguish QED corrections from higher dimensional operators: both contri¬ 
butions give rise to higher moments but crucially the QED corrections dominate for 
moments of increasing Q, cf. the discussion at the end of section 5.3 and specifically 
equation (65). A J 7 -wave at the amplitude level contributes to a Q = J 7 + 1 moment 
through interference with the SM Pg- wave. We conclude that QED and higher spin op¬ 
erators could be comparable for but for G^ >3 - one would expect the former to 
dominate. 16 

14 A partial result, photon emission from initial and final state, was reported in [60]. 

15 Collinear photon emission in the inclusive case was studied recently in [61]. The additional photon of course 
leads to terms which go beyond the Pff angular distribution. Note, in view of the presence of these terms 
through virtual corrections they also have to be present in real emission by virtue of the Bloch-Nordsieck 
QED infrared cancellation theorem [62]. The authors [61] find within their approximation that the third and 
fourth moment are two orders of magnitude smaller than the leading contributions. This is in the expected 
parametric range but one cannot draw precise conclusions on the size of this effect for the exclusive channels 
discussed in this paper. 

16 Another criterion could be that corrections from higher spin operators are uniform in the lepton mass provided 
that lepton flavour universality is unbroken. This is though delicate since the measurement of Rk questions 
this aspect. 
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• Lepton-flavour dependence of QED corrections: differences between and Ggf~ 3 ^ 

in the range above q 2 > 1 GeV 2 indicate the importance of the flavour dependence. This 
gives an indication on how much the branching fractions (zeroth moments) and therefore 
Rk is affected by QED through lepton mass effects. Note that due to 1 r 1 ( to ^ / m {,) - e ffec ts 
it is conceivable that is small, say 0(1%), but that Oge~'^ is larger. Note for 

example that G/j]l = Afb is consistent with the SM prediction excluding QED, which is 
O(m^), within errors in the few percent range [37]. 

5.4.2. Combinatorial background in B —» K* /i+ fi~ below the narrow 
charmonium resonance region 

A characteristic feature of B —> K^£ + £~ transitions is the large contribution to the branching 
fraction through the intermediate narrow charmonium states J/T and T(25). For example 
B(B KJ/'fr)B(J/'S> -> p + p~) — (8 • 10 _4 )(6 • ICE 2 ) ~ 5 • ICC 5 is three orders of mag¬ 
nitude larger than the measured differential branching fraction, dB(B —> K*/i + p~)/dq 2 ~ 
2 • ICC 8 / GeV 2 [63], well below the narrow charmonium resonances region. It is therefore 
legitimate to be concerned with possible combinatorial backgrounds in this region. 

Assuming that such backgrounds are relevant this raises the question as to how they can 
be distinguished from the signal event. In the case where they can be absorbed into the 
background fit-function they would not impact on the analysis. Whether or not this is the 
case is a non-trivial question. Pragmatically, however, background events can be expected 
to perturb the hierarchy of the moments as compared to the true signal event. One would 
expect the background events to fall off only slowly for higher moments in the lepton partial 
wave. 1 ' Hence the size of these effects can be diagnosed through the measurement of higher 
moments as a function of q 2 , independent of model assumptions. By the latter we mean that 
higher moments peaking below the charmonium resonances will be indicative of the type of 
combinatorial background mentioned above. 

A possible example of such backgrounds is the process B —> Kpfl 7 where the photon is 
not detected but energetic enough to cause a significant downward shift in q 2 = (£\ + tfl) 2 ■ 
Such an event would be rejected as a B —> Kpflp,~ signal because the reconstructed B-mass 
tur^ would fall outside the signal window (i.e. vuk^ < m# — A and A ~ 0(100 MeV)). 
If additionally a 7 r-meson from the underlying event is detected, the event could conspire to 
enter the signal window of B —> K*{—> Kir)p, + p,~ (i.e. tur^^^ — mg and toav — tur*). It is 
therefore conceivable that the small chance of the events described above is overcome by the 
enhancement by three orders of magnitude of the J/T transition. If such events are present and 
not rejected then this leads to a bias in B —> K*{—t Kn)p, + fi~ transitions below the narrow 
charmonium resonances. More precisely, denoting the momentum of the undetected photon 
by r, the shift in q 2 is as follows q 2 ~ m 2 = (^1 + £2 + r ) 2 <? 2 ignal = (h + 12) 2 < ™ 2 j/y 

This is particularly relevant as some of the anomalies from the LHCb measurements, in 
particular the angular observable P 5 , are most pronounced in bins just below the J/T-resonance 
[26,27]. To what extent such operators correspond to new physics in O 9 = Oy [64,65] or effects 
from charm resonances [ 66 ] is a difficult question since they contribute to the same helicity 
amplitude. They can be distinguished from each other by analysing the g 2 -spectrum of the 


1( Similar things can be said about the hadronic partial wave, but as the detection of the Pk -wave is part of 
the signal selection the presence of such higher waves would have less influence. However, the remaining 
background might impact on the Sk- wave, which does matter since the Sk -wave enters the analysis. 
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observables and by the determination of the strong phases which can originate from the charm 
resonances [66] . This could be through the determination of the complex-valued residues of the 
resonance poles [66], or simply the strong phase in the region below the g 2 -resonance through 
Im[Gi ’ ] ~ Pq, which corresponds to the imaginary part of Re[G^ - 1 ] ~ P,i (33). 


6. Conclusions 

In this work we have generalised the standard helicity formalism to effective field theories 
of the b —> s££- type. The framework applies to any semi-leptonic and radiative decay. The 
formalism has been used to derive the angular distributions I^l (28) and I® (39) for non¬ 
equal lepton masses with the full dimension-six effective Hamiltonian, including in particular 
scalar and tensor operators. Explicit results for B — > K*£\£^ and for B —> Kt\l^ can be found 
in appendices C and D respectively as well as a Mathematica notebook (notebookGHZ.nb) 
provided in the arXiv version. Comments on differences conversion of observables between 
theory and experiment with the literature are reported in appendix C.2.1. Minor discrepancies 
in tensor contributions with respect to previous results are discussed in appendix C.1.2. 

The approach clarifies how the lepton factorisation approximation determines the specific 
form of the angular distributions I^l and I^\ and how these distributions are extended by 
the inclusion of virtual and real QED corrections, as well as higher-spin operators in the 
effective Hamiltonian. Higher-dimensional spin operators provide new opportunities to search 
for physics beyond the SM. We have argued that, within the SM, QED effects and higher- 
spin operators can be distinguished from each other by their differing fall-off behaviour in 
increasingly higher moments in the f^-angle. 18 

Assessing higher moments can shed light on current anomalies with respect to the SM. We 
have argued (cf. section 5.4.1) that higher moments in B —> K£ + £~ {£ = e, //) are a window 
into QED corrections and therefore of importance with regard to the Rk measurement [25]. 
In view of tensions of angular predictions in B —> with experiment [26,27], the 

higher moments can be of help in assessing their origin, such as the possible leakage of J/T 
events into the lower nearby q 2 -bins (cf. section 5.4.2). As another example we mention the 
R(D^) = B(B -> £>(*) tu)/B(B —> fiv) ratio measurement [67-69], suggestive of some 
tension with the SM. A higher moment analysis could again be useful in assessing the impact 
of QED, lepton mass or cross channel backgrounds on these results. 

To measure and bound higher moments is relevant as their contributions can bias likelihood 
fits. We therefore encourage the investigation of higher moments in several experimental 
channels from the various perspectives discussed above. 
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18 In addition higher-spin operators can be distinguished from QED corrections by universality in the lepton 
flavour. However, it should be kept in mind that lepton-universality is questioned by the Rk measurement. 
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Note added While this paper was in its final phase a paper using the helicity formalism for B —> 
KH+f.- appeared [70]. The paper uses the standard Jacob-Wick formalism and therefore includes HAs 
of definite spin. This is an approximation that holds up to lepton mass corrections in the SM and does 
not allow the inclusion of scalar operators for example. 

Changes in conventions and presentation Notational changes with respect to the arXiv 
version 1 , aimed at clarifying the underlying structure, are as follows: i) results are presented for 
B — > rather than the conjugate decay, ii) we use C^ Wilson coefficients in place of Ct( 5 ) for 

the tensor operators cf. appendix A.2 for details, iii) the angular distribution (C.l) is presented in terms 
of gt in place of Jy in order to emphasise the differences of angular convention of this paper and the 
theory community (as discussed in appendix C.2), iv) lepton HAs are presented in the A , V rather than 
L, R basis and v) timelike HAs are absorbed into scalar and pseudoscalar HAs. In addition we provide 
a Mathematica notebook, entitled notebookGHZ.nb, containing the results presented in appendix C.4 
for the decay mode B —t K*(—t Ktt)£iI 2 for non-equal lepton masses. 


A. Results relevant for all decay modes 


A.l. Decomposition of SO( 3,1) into SO( 3) up to spin 2 

The aim of this appendix is to give some more detail about the decomposition (15) and in 
particular extend it to the two-index case, which includes the discussion of spin 2,1, 0. 

In section 2.2 it was shown that insertion of the completeness relation (15) corresponds to 
the decomposition, or branching rule, 


(1/2) 1/2)50(3, 1} 


50(3) _>(1 +3)50(3)5 


(A.l) 


where (1/2,1/2) is the irreducible vector Lorentz representation. We remind the reader that 
the irreducible Lorentz representations, denoted by (./1 ,. 72 ), are characterised by the eigenvalues 
of the two Casimir operators of SO( 3,1). Inserting the completeness relation twice therefore 
corresponds to taking the tensor product (1/2,1/2) 0 (1/2,1/2) which decomposes as 


SO(3) 


((1/2,1/2) ® (1/2, 1/2)) so(3i1) = ([(1,1)] © [(1,0) 0 (0,1)] 0 (0,0))so ( 3,i) 

vs ytj j 

([1 • 5 0 1 • 3 0 1 • 1] 0 [2 • 3] 0 1 • l) so(3) = (l-503-3©2 - l) SO (3) ■ (A.2) 

The double completeness relation 

da/39^8 = ^ap'yS + SafaS + ^a/3jS (A.3) 

can be decomposed 


2 j 


^a/3^S 




'y ] y ^ w cry w /3 8 ’ u a0'y8 


1 1 

A \ ' , ^tt 




E , rtt _ , ,tt —tt 

> Oa0"/8 ^a'y^pS ) 


J=0 X=-J 


A=—1 


A=—1 


into parts containing zero, one and two timelike polarisation vectors 


(A-4) 


= cj Q (t)o; 7 (A) , = o; Q (t)o; 7 (t) , 
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W «7 — X] C , AAi 1 A2 Wq ('^ 1 ) U; 7(^2) , 
Ai, A2 = 1 


(A.5) 


with A = Ai + A2 in the first term and the polarisation vectors Lo a ( A) are parametrised as 
w M (±) = (0,±l,i,0)/V2 , ^ M (0) = (q z , 0, 0, go)/ \[q* , u/*(i) = (go, 0, 0, q z )/\fq 1 , (A.6) 


which we reproduce from (16) for the reader’s convenience. A few explanations seem in order. 
The minus sign in front of g in (A.4) is due to there being an odd number of timelike 
polarisation vectors. The first, second and third term in (A.3) correspond respectively to 
the (1,1)-, [(1,0) © (0,1)]- and (0,0)-terms in (A.2). It is convenient to rewrite the double 
completeness relation (A.3) in a form that makes the decomposition into the different spins j 
explicit 

2 j 

9a.p9y8 = e a7 ■ e f35 ■ (A-7) 

J=0 X=-J 

Above the scalar product stands for 


A -j, A' _ r 
^cz~/'^/35 *J0 


. , 0 , 0 pr 0,0 | tt -tt 


(3S 


+<5ji 


1 A—1,A' t,\ —t,A' t.X — t,X' 

^a-y^pS W 7Q: W <5/3 


+ <5j2 


,2^2 ,a' 

J Q7 U P8 


(A.8) 


The single completeness relation (15) in the analogous notation of (A.7) reads 


1 


J 


9ap — / 4 

J=0X=—J 


J,A ? J,A 
/ , e a e p 


(A.9) 


with e J a X = 5j±uj a (\) + 5j 0 uj a (t). 

When applying the double completeness relation to generic decay structures, it can be seen 
from (A.4) that in general one expects two distinct contributions to the amplitude from S t a j 3ryS , 

H IIV L' W -> - (UtxCtx + UxtCxt) + ■ • • , 

where Jitx = and analogous notation for 7i\t, Ctx and C\t- If, however, the objects 

H^ u and U J,V are both symmetric or antisymmetric in the Lorentz indices, then T-L\ t C\t = 
l~Ltx£tx and the two contributions can be combined. We have used this simplification in defining 
the generalised HAs for the B —> K*£ ±£2 (23) and B —> K £\£.2 (36) decays respectively, 
resulting in the extra factor of 2 associated with the terms hJ' C7^ ^ h Tt C^ relative to 
other contributions in the generalised HAs. 

A.2. Additional Remarks on effective Hamiltonian 

Here we collect a few additional remarks to the effective b — > s££ Hamiltonian quoted in 
Eqs. (12,13). Contributions proportional to have been neglected. The chromoelectric 

and chromomagnetic operators O7 and Os, along with the contributions of the four-quark 
operators Oi,...,6) can be absorbed into Oy through defining an effective Wilson coefficient 
= Cg ff . We can rewrite = 1/2 (Ot ± Ot 5 ), with the latter defined as 

O t = sa^ u b£a^ v £ , O t 5 = sa^^bla^ , (A. 10) 
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(note: Ot 5 = sa^bia^'y^i = — |e Q/ 3 /il/ sa a pb£a pu £ with the last equality depending on conven¬ 
tions) and the relation between the Wilson coefficients is therefore 


Cj = Ct± Ct 5 , Ct( 5 ) = -(CV ± C ' T ) (A. 11 ) 

in the sense that CtOt + CtsOts = Cq-Oq- + C' t O' t . 


A.3. Definitions and Results of Leptonic Helicity Amplitudes 

The calculation of the Leptonic Helicity amplitudes is an important part of the generalised 
helicity formalism described in this paper, and the method for their calculation is outlined 
in [3]. In the Dirac basis of the Clifford algebra, with ex* as the usual 2x2 Pauli matrices, 



the particle u and anti-particle v spinor are given by 

= (V^i +me i, 0 , \JE\ - 7714,0) = {fit’ 0 ,/ 3 r, 0 ) T , 

= (0, \JE\ +m ei ,0,-y/Ei = (0, , 0, -/3f) r , 

v = (y/E 2 -mi 2 ,0,-y/E 2 + m e2 ,o') = (/FT, 0, 0) T , 

v (“0 = (°, y/ E 2 ~ me 2 , 0 , \Je 2 + my) = ( 0 ,^, 0 ,/ 3 ^) T , 


with implicit definition of /3- = y / £'i ± my. The spinors are normalised as ti(Ai)w(A 2 ) = 
<!>AiA 2 2my and u(Ai)?;(A 2 ) = — b\ x \ 2 2mi 2 . The leptonic HAs (18) contracted with polarisation 
vectors give rise to the HAs £aiA 2 


££a 2 = (^i(Ai)f2(A 2 )K T A >|0) = h(Ai)T x u(A 2 ) , (A.13) 


(where t\ = e~ for example) and the F x |a x ->.a^ (A^ = Ai — A 2 ) as defined in Tab. 3.1. Using all 
the equations above the evaluation of the lepton HAs is then straightforward and the results are 
presented below, for lepton masses my 7 ^ my in the first set of matrices and my = my = mi 
in the second set . 19 The first row (column) corresponds to Ai(A 2 ) = — ^ and the second row 
(column) corresponds to Ai(A 2 ) = +|. For the B —» K*l if 2 decay mode, ie t\ = £~, the lepton 
HAs are given by 


£ y (Ai,A 2 ) 



-y/2 (/?+# +/?r/? 2 -) \ 

( 2me 


V2 (/3+P+ + /3f/3 2 -) 

PfPt - f^fr 2 ) 

^ V -v'j? 

2 mi J 


19 The expressions for mt 1 ^ me 2 can be applied to studies of lepton flavour-violating processes in all the decay 
modes considered in this paper within the lepton factorisation approximation, and are also applicable to 
decays involving an W in the final state e.g. B —» D*lv. 
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0 yfi&Pl 

~VW 2 Pe o 

0 

0 \fq 2 Pi 

y/q 2 0 

0 -y/q* 

-iy/WPt 0 

0 i^/UpPg 

iyfq 2, —2iy/2mg 

—2iy/2mg iy/f 
(A.14) 

where pf 2 = ± I rng 1 2 as before. Above we have used (3~ —> E/3i for 2 -> since 

FI 2 = g 2 /4, where E is the energy of either lepton in the rest frame of the lepton pair. Note 
that the scalar transitions S and P are necessarily diagonal since Xg = X\ — X 2 = 0. Timelike 
vector and axial Lepton HAs are integrated into the Hadron HAs (C.15). 


Pi Pi ~ Pi Pi \/2 {PtP 2 + Pi Pi) 

-V2 {Ptp 2 + pi Pi) pi Pi - Pi pi 

Pi Pi + PI Pi 0 

0 Pi pi+pi Pi 

PiPi +pi pi 0 
0 -ptPi -pipi 
- w 2 {ptPi +pi Pi) -2 i {pfpi - pi pt) 

2% {Pt pi - pi pi) w 2 {pt pi + pi p+) 

i {pi pi + pi pi) -w 2 {pi pi - pi pi) 

-iV 2 {pipi - PiPi) i {pipi + PiPi) 


B. Details on Kinematics for Decay Modes 


While within the formalism described in this paper it is not essential to consider the full 
kinematics of the decay, as the evaluation of the hadronic and leptonic HAs can be performed 
within their respective rest frames, we collect here the kinematics used in calculating the 
angular distribution using the Dirac trace technology approach [22, 23] in order to facilitate 
comparison. The Kallen function that often appears in our results is given by 

A (a, b, c) = a 2 + b 2 + c 2 — 2 (ab + ac + be) . (B.l) 


For a decay A —> B + C, in the rest-frame of A, it is related to the absolute value the spatial 
momentum of the B and C particles as 


\Pb\ = \pc\ 



m c) 


2 mA 


(B.2) 


B.l. Basis-dependent kinematics for B —> K*£i£ 2 

We parametrise the kinematics of the (£\ = £~ and £2 = £~) 

B K *{->■ K(p K ) 7 r(p w )) h{h)Hh) (B. 3 ) 

decay mode. To do so we need all four momenta p n , px (p = Pn + Pk), £1 and £2 {q = £1 + £2) 
in a specific frame for which we choose the H-restframe. It is simplest to first obtain £1.2 and 
p n ,K in the restframe of the lepton pair and the A'*-meson respectively: 

^1,2-restframe : £ 1 = (Ei, \p £ \ £) , £% = (E 2 , - \pg\ £) , 
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and the definitions 


where 


— \ f„ 2 


= (E k , \pk\ k) , 

Pn = (E-k, ~ \pk\ k) , 

(B.4) 

sin</>sin#£,cosf^) , 

“ ■ 5 ■ 


cos 9 k) , 

1^1 V Aa'* 

\Pk\ = 

ZTYIk* 

(B.5) 

X(m 2 K *,m 2 K ,ml) , 

X B = X(m 2 B ,m 2 K *,q 2 ) , 

(B.6) 


are the explicit Kallen function used throughout. The lepton and hadron energies are then given 

by E \2 = \J m J 12 + \Pe\ 2 , En,K = + \pk\ 2 and obey E 1 + E 2 = \fqt 2 and E n + E K = 

m K * ■ 

The polarisation vectors t/^A) of the Ii*-meson in its restframe, using the convention in [18], 
20 


are 


(B.7) 


rf(0) = (0,0,0, 1) , rf(±) = (0, Tl, i, 0)/V2 . 

In the B-restframe, p B = (m#, 0,0,0), the momenta take the following form 

(hY = {fi(Ei, q 0 , q z ), \p e \ sin 6 e cos 4>,-\pi\sin 9 e sin 4>,f £ (Ei,q z ,q 0 )) , 

(hY = (fe(E 2 ,qo,-qz),- |^| sin cos (/>, + \pi\ sin 9 e sin fe(E- 2 , q z ,-q 0 )) , 

(. PkY= (fK*(EK,Po,qz),-\PK\sin9 K ,0,-fK*{E K ,qz,Po)) , 

{p-kY = ( fK*( E ir>Po,-qz),\PK\sin9 K ,0,-f K *(E 7T ,q z ,-po )) , (B.8) 

with fi{a,b,c) = {ab + c\pi\ cos 9 1 )/and /at*(o, 6, c) = (ab + c\pk\ cos 9 k)/ m* K , and it is 
easily verified that 

q^ = (4 + hY = (qo,0,0,q z ) , P^ = {pk +PnY = (po,0,0,-q z ) , (B.9) 

(po = Ek*) while the polarisation vectors of the K* in the B-restframe are 

<(0) = (-q z ,0,0, Po )/m K * , ^(±) = (0,^l,i,0)/v / 2, (B.10) 

where Po + qo = m B and q z = \fXs/(2m B ), in accordance with (B.2), is the three-momentum 
of the lepton pair. 

For completeness let us add that in the case of: 

• B —> K*l\{li)l 2 (l 2 ) the replacement rule i —> i </,_>_</, = (cos^sinf^, + sin</>sin^, cos9i) 
applies. Note this is coherent with figure C.2 in the next section; 

• identical lepton masses the following replacements are in order 

Ei, 2 V^/2 , t/V -> (<? 2 )& (B.ll) 


/ AtY) ^ 

where we recall that = y 1- ■ 

20 The polarisation vector r\ corresponds to 7 in [18] (c.f. appendix A therein). The exact correspondence between 
the convention used in [18], and also in this paper, and the Jacob-Wick convention [1,3] is 7 (±) M | j-^gj = 

—77(±) m I |ij, 7?(0) m | | 7 gj = 7 ( 0 )^ jij. The final distributions remain the same but the off-diagonal elements 
of the lepton HAs (or matrices) change sign (A.14). Note in particular that the hadron HAs (C.15) remain 
unchanged. 
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B.2. Basis-independent kinematics for B —> K*£ 1 £ 2 

Introducing the notation 


q» = - i 2 y , p» = {p K - Ih y, 


(B.12) 


in addition to (B.9). the invariants that can be formed out of p, P, q and Q are given by 


q ' Q = ml - rnj 2 , Q 2 = 2 [m\ + mj 2 ) - q 2 , q ■ p = - (rn 2 B - m 2 K * - q 2 ) , 


p. P = m 2 K — m 2 , P 2 = 2 (m 2 K + m 2 ) — m 2 K * , q - P = 


2p ■ P q ■ p + cos OkV^B^K* 
2m 2 K „ 


p ■ P y / AsA 7 * cos 9g + 2q ■ p^/\k*\* cos 6k cos 9g yf A k* A 7 * sin 6k sin 6g cos 4> 

Q ' P = n — 2 2 


2m z K «q 


m K - 


V? 


+ 


Q-p = 


q - Qq - P 


2q ■ Q q ■ p + cos 6g y / AsA 7 


t(P,p,Q,q) = - 


2q 2 

sin sin da sin <^^/AbAx*A 7 

2 m K *\fq 2 


(B.13) 


with p 2 = m 2 ^*, the A’s defined in equation (B. 6 ), e(P,p,Q,q ) = e a p 1 sP a p l3 Q' y q s and the 
6 oi 23 = 1 convention for the Levi-Civita tensor. Note that the kinematic invariants for B —> 
K *are the same up to e ( P,p , Q , q) —> —e (P,p, Q, q) which originates from the only 
change in angles 0 —> — cj). 


C. Specific Results for B —> K*(—y K7r)£i£ 2 


C.l. Fourfold Differential Decay Rate 

The angular distribution for B —> K*(—t Kk )£\£2 is usually presented in the form (e.g. [35]) 


8 vr d 4 T 

3 dq 2 dcosdi dcos&K dcj) 


j(0) 

-yp = (gis + 92 s cos 26 g + g 6s cos 6g) sin 2 d K + 

(gic + 92c COS 26 i + gee cos 9?) cos 2 d K + 
(53 cos 2 4 > + 59 sin 2 (f>) sin 2 6k sin 2 dp + 
(54 cos 4> + gs sin <£>) sin 26k sin 29£ + 

(5,5 cos 4> + 97 sin 4>) sin 26k sin 9g , 


(C.l) 


which can be condensed as 

~tt T 77 - 7~7 — n—72 = Re [ ( 91s + 92s cos 26g + g 6s cos 6p) sin 2 d K + 

3 dq z dcos6g acosP/6 dep 

(,gic + 92c COS 29£ + 56c COS 6g) cos 2 6 k + 

g-2*</>g 3 gj n 2 q k g - n 2 _|_ 

sin 26k {Qi sin 2 6g + Qg, sin 6g ) ] , (C.2) 
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where we have defined 


£ 3 , 4,5 = (53,4,5 + * 59 , 8 , 7 ) • (C.3) 

We have introduced the notation gi rather than J* in order to minimise the potential of confu¬ 
sion due to the angular conventions discussed in appendix C.2. The relationship between the 
gi(q 2 ) and the Gm’ 1 * ( q 2 ) was given in (32) but is repeated here for convenience: 

Gq '° = g (3 (5lc + Igis) ~ (52c + 2$2 s)) , Go’ 1 = - (56c + 2 c/6s) , Gq’ 2 = — (52c + 2<72 s) , 

Gq ° = ^ (6 (5ic - gis) - 2 ( 52 c - 52s)) , Go" 1 = ^ (g 6c - g 6s ) , G 2 0 ’ 2 = y ( 52 c - 52s) , 

G 2 ’ = - 7 = (55 + * 57 ) , G 2 ’" = — (54 + * 58 ) , G 2 ’ = — (53 + * 59 ) • 

V 3 v -v-" 3 v- v -' a v- v -" 

=05 =04 =03 

(C.4) 

Explicit results for the Gm’ le are presented in section C.3 for the case of identical final-state 
leptons m£ 1 = me 2 and section C.4 for the more general case m / mi 2 . 

C.1.1. Kinematic endpoint relations in terms of G l ^ ,le 

In Ref. [18] it was shown that the HAs obey symmetry relations at the kinematic endpoint 
due to symmetry enhancement. This is due to the K* being at rest resulting in symmetry in 
all space directions i.e. helicity directions. The relations for the HAs in equation (13) in [18] 
lead to the following equivalent of equation (21) in [18] 

Go° + 0 , G 2 ’ 2 Re[G 2 ’ 2 ] , G 2 ’ 2 -> -2Re[G 2 ’ 2 ] , G 2 ’ 2 2Re[G 2 ’ 2 ] , (C.5) 

with all other five Gm’ le vanishing. Recall that Gq 0 is proportional to the total decay rate. The 
relations between the GJi are not accidental but have to do with the symmetries of a multiplet. 
The factor of two between Gq' 2 and G 2,2 2 y once more, originates from absorbing G 2,2 ^ into 

G 2 ’ 2 ). The results of the threshold expansion, linear in the K* momentum k ~ X(m 2 B ,m 2 K *, q 2 ), 
can be inferred from equation (30) in [18] taking into account the different angular conventions 
detailed in figure C.2. 

C.1.2. Comparison of angular distribution with the literature 

The angular distribution (C.l) was first computed in the SM for massless leptons in [22], 
extended to include equal lepton masses in [23,71]. A full dimension-six operator basis was 
considered in [72]. The basis was extended to lepton mass corrections for (pseudo) scalar 
operators in [34], enforcing the g6c-structure, and tensor operators by the authors in [35,73]. 
We compare our results with regard to [35], which is the latest reference. 

Taking into account the change 54 , 6 , 7,9 —> ~ £ 4 ,6, 7,9 (cf. figure C.2) and comparing at the 
level of form factors (naive factorisation) only we find agreement except for tensor interference 
terms. Agreement is established when Ct 5 —> —Ct 5 in [35]. The latter might be related to 
the fact that the relations tr^q^y'^ys] = 4A ie 0 ^ 5 and fJ a/3 75 = ~X^e al3 ' y5 a 1 s (with A = ±1 
depending on conventions - A = 1 in this paper) are not consistent with equation C.16 [35] 
(v3). 
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A minor difference is that the authors of [35] have chosen not to present the tensor contri¬ 
bution in Js,9(<78,9)> since such contributions vanish in the narrow-width approximation. 21 In 
addition, we find that a few of the HAs in [35] equation C.13 do not agree with their definitions. 
These disagreements do, however, drop out in the final expression. 


C.2. Angular conventions 


In this section we discuss and compare the LHCb and theory angular conventions. The main 
result is shown in form of a commutative diagram in figure C.2. We proceed by first discussing 
the CP-conjugate modes in each case and then link the conventions with each other. 

The LHCb conventions [31], which are the same as adapted in this paper, are shown in figure 
1. The rationale behind the definition of the conjugate mode is as follows. Firstly, particles 
are mapped into antiparticles, corresponding to a C-transformation. Then the momenta of all 
particles are reversed, changing the angle cf> —> 2ir — </>. This leads to sign changes in 57 , 8 , 9 . 
Hence the conjugate mode corresponds to a full CP-transformation 


d A r 


d 4 r 


dq 2 d cos 6id cos 6Kd(/) LHCb dq 2 d cos 9id cos 9xdcj) 


CP 


and the quantity 


d 4 (r±f) 


dq 2 d cos 9id cos 9xdcj) LHCb ’ 

is therefore CP-even (-odd). Above T = F(B — > K*t \l-i) and f = F(B —> 

The theory conventions for CP conjugates are such that they facilitate the implementation 
of decays which are not self-tagging (such as B S ,B S —> </(—> K + at hadron collid¬ 
ers). When going between conjugate modes the conventions are that the angles transform as 
(9i,9x,(j)) — > (-7T — 9i, it — 9k, 2vr — 4>), 22 which leads to sign changes in 55 , 6 , 8 , 9 - This transfor¬ 
mation rule corresponds to a full CP-conjugation, but with the angles 9i, 9 k associated to the 
same particle rather than the anti-particle. 

To find the transformation between the theory and LHCb conventions is not straightforward 
because it is difficult to find a theory paper that resolves the four-fold ambiguity of defining the 
angle 4> and/or shows a figure consistent with the definitions used in the corresponding work. 
We have taken a pragmatic route in verifying that the results in [34,35, 71] agree with each 
other for common contributions, and crucially that our results are in agreement with these 
contributions for B —> K*t\t2 if J4, 6 , 7,9 = — 54 , 6 , 7,9 and Ji, 2 , 3 , 5,8 = 51 , 2 , 3 , 5 , 8 - This completes the 
diagram in figure C.2. 23 

In summary, 

d 4 { F ± f) 

LHCb 


dq 2 d cos 9id cos 9 xdcj) 


4A 


S[A] 1 , 2 , 3 , 4 ,5,6,7,8,9 , 


21 In v3 of [35] it is stated that agreement with v4 of [74] is found up to a sign of an interference term between 
a scalar and a tensor HA. This suggests that we agree with [35] but disagree with [74] on that sign, as well 
as the sign of Cts- 

22 Equivalently one can use the angular redefinitions (n — 9 (, 9 k , 7 r — (j>) and ( 9t — Tv,9K,2n — (ft), which are 
sometimes stated in the literature. 

“ J One can come to the same conclusion in another way [75]. Let us again consider B —> K* £ 1^2 ■ In general 
9 k I LHCb = 9 k I theory is chosen to be the same angle and the theoretical community chooses 0r|LHCb = 
7T — 9 1 1 theory ■ The only unknown remains the angle cf>, for which one may use the scalar- and cross-product 
definitions. Using Appendix A in [31] and likewise in [71], we infer that cos0|LHCb = cos (^theory and 
sin0|LHCb = — sin ^| theory- Taking all angular changes into account this results in sign changes in 54,6,7,9 
which is consistent with our explicit computations mentioned above. 
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QA 5 9 —^ —J 4 5 Q 

B^K'U) LHCb —^B-► tf*^| theory 


#7,8,9 —»• —#7,8,9 


^5,6,8,9 —>• —^5,6,8,£ 


B -> K*^| L HCb 


#4,6,7,9 —>• —^4,6,7,£ 


5 —> if* ^(theory 


Figure 4: Changes of angular functions when going from one mode to the other. For CP conjugates 
the conjugation of the CP-odd (weak) phases are suppressed. Angular functions whose signs 
do not change are not indicated. 


d 4 (r ± r) 


o *S'[ J 4]l,2,3,4,7) ^4[<S']5,6,8,9 , 


dq 2 d COS 9id COS 9xdcj) theory 

where the CP-even (-odd) quantities are 

Si[Ai] = y + r CP 

with adapted notation from [34], Written in yet another way (C.6) is equivalent to 


9i ± 9 


CP 


(C.6) 


(C.7) 


(flS A*5)l,2,3,5,8|LHCb = + (J, A, S)i t 2,3,5,8|theory , 

(g, A, 5')4 i 6,7,9|LHCb = “ (#, A *5)4,6,7,9|theory • (C.8) 

In order to understand (C.6) and (C.8) one has to keep in mind that B -A- K*{1\I 2 ) rather than 
its conjugate is the reference decay. Note that at the LHCb (hadron collider) B s —> cj)g + g,~ is 
untagged and therefore, setting aside the issue of production asymmetry, only <S'i, 2 , 3 , 4 , 7 5 A 5 , 6,8,9 
are experimentally accessible. 


C.2.1. Angular observables in the literature and conventions 

//\ 

We aim to find the relation between angular P ?; ' observables as defined by the theorists [33] and 
their adaptation by LHCb [26] . In matching the results and creating the dictionary one needs 
to pay attention to the fact that [26] and [33] define the P\ in terms of g% and J,: differently, as 
well as the different angular conventions for gi and per se (shown in figure C.2) ! 

Amongst the twelve observables discussed in section 3.2, eight of them, Pi, 2 , 3 , P 4 5 6 8 an( 4 
Afb> depend on angles and definitions. 

The P[ and Afb are defined by LHCb [26] as 


-^4,5,6,81 LHCb 


*54,5,7,8 [LHCb 


71fB |LHCb = 


3 (Sq s I LHCb) 

4(r + f) 


(C.9) 
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where S) is defined in equation (C.7) and 2A/^ in = \/Fi( 1 — F in our notation used in section 
3.2. LHCb has not defined Pi 2 3 and we shall assume the same functional form as in the theory 
paper [33]. 

In [33] the eight equivalent angular observables are defined as follows 24 


Pi 


P2 


Ps 


p' 

-*4 

P' 

-'5 

Afb 

P' 

-'6 


p' 

r 8 


1 


2S 2s 

1 

8 S 2 S 

-1 

45^ 

1 

•^bin 


5 3 


S es 


S 9 


s 4 




S 5 


3Sq s 


4(r + r) 




-1 

■^bin 


^8 


TTTi ('“'3 | LHCb) — +-Pl|LHCb> 

^gT“(-'S'6s|LHCb) = |LHCb i 


4 g 2 ( S 9 1 LHCb) 

= — P 3 LHCb j 

K n ( S 4 LHCb) 
•'Win 

= — 2P4| L HCb 

9A r/ (^l LHCb) 
Z - /V bin 

= + P 5 1 LHCb , 

3( — Sq s LHCb) 

= +A FB LHCb 

4(r + f) 

9A r/ ( 571 L HCb) 

zyv bin 

= + Pq LHCb , 

K r, (^8 1 LHCb) 
- /V bin 

= - 2Pg|LHCb 


where we have directly translated into the LHCb conventions. It seems that we differ from the 
theory community in the sign of the observables S’[^4] 7 , 8 , 9 - For example, both Pg = Pg|LHCb 
and Pg = — 2Pg|LHCb differ from the relations given in the caption of table 1 in [64] by the 
aforementioned sign. Our relation iS[A]g = — 5'[H]g|LHCb also differs from the one given by [65] 
in table 1 by a sign. 


C.3. G l ^ tl for B —> K*£ 1 £ 2 in terms of Helicity Amplitudes for m<. = me 

When the masses of the leptons are identical, we obtain for Gm’ le = Nq 2 G l £' h (with A f defined 
in (26)), 


<° = | (1 - m/) (| Hi | 2 + |PlT + |P 0 T + (V -> A 


rV,2 


rV,2 


+ 


4m/ 


Hl \ 2 + |P y | 2 + |PT^| 2 - (V -> A)) + |/31 |P S | 2 + | |P p | 2 


+ l (1 + 8r «^) ( 


H 


T t 


+ 


H 1 


+ 


H Tt 


n + ^[\hI\ 2 + \h t _\ 2 + \h p \ 2 


16 * T 

+ —me Im 

O 


Gg’ 1 = (Re [P+P+ - H V H^] + Im \V2H%H p + 2 H Pt H s 


P+P+ + H v H Pt + H^H Pt 


— 2/rqRe [Hq H s ] + 4m^Im 


Hi H r < - H*H 


AfrTt 


Note that [34,35] define ^fb = 4 ^r+r) which results in Afb| [34 35 ] = —Hfb|LHC b- 


[34,35] 
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< 2 = 


2 H^\ 2 - \HX\ 2 - Hi * + (V -> A) - 2 (2 H'i ‘ * - Hi 


T v i 2 


r T|2 luTI 2 


r T|2 


-4 2 




Hi 


Hi 


Gq’° = (! - m* 2 ) (|^| 2 + \H V _\ 2 - 2 K| 2 + (V A)) 


(Kl 2 + \Hi\ 2 - 2 Kl 2 — (V —> A)) + |/3| |F 5 | 2 + ^ |F P | 2 


-g(l+ 8 ™e 2 ) ( 


H 


T t 


+ 


H 


T t 


- 2 


H Tt 

n o 


2 \ 4 

9 ; 


j2 ^ I Hi 1 2 + I Hi 1 2 — 2 I 1 2 


16 ~ T 
- 7^ m t -*- m 


Gq 1 = (Re [HXHX - Hi Hi] - 21m \V2 H^ H p + 2 H^ H s 


W + HX + Hi Hi! - 2 Hi H ( y 


+ 4 mi Re [Hq H s ] + Im 


/■/;' h 'X - Hi Hi 


G o 2 = - 2 ( /l ( 4 K| 2 + \K | 2 + \Hi\ 2 + (V^A)~ 2(4 K| 2 + |^| 2 + \Hi \ 2 


-4 4 


H o 


+ 


H? 


+ 


H 


Tt 


G 2 ’ 1 = ^ (HXHi + H^hX - - Hi Hi + 2m e ( HXH S + H s Hi) 

V3 


- V2i (H p Hi - H p H p + V2 (H s Hi - HXH S 


- 4irh e ( H^hX + H^Hi - H^H^ - H^Hi 


G 2 ’ 2 = ( Hi Hi + Hi Hi + (V -> A) - 2 ( Hi Hi + Hi Hi + 2 ( H'^HX + H^Hi* 

O 


fTt rjTt _L_ U T t TjTt 


G 2 ’ 2 = [Hi Hi + Hi Hi - 2 [Hi Hi + IH^Hi 


tV fjV 1 rr.4 fVA 


t T fjT 


T T t rjTt 


(C.10) 


where and we recall that /3^ = y^l — 4m 2 . The number m in Gm’ l( corresponds 

to the units of plus helicities. The common factor of g 2 in all observables as compared with 
standard literature results is a consequence of our choice of normalisation, whereby all global 
factors are placed outside the HAs. The factors of i where they appear (explicitly and implic¬ 
itly) in G 2 ’ 1 , G 2 ’ 2 and G 2 ’ 2 are not accidental, as the results given above are complex and one 
must take the real and imaginary parts of these results to recover the observables 53 , 4 , 5 , 7 , 8 , 9 - 
Note that it is sometimes convenient to express results in terms of the transversity ampli¬ 
tudes, which possess a definite parity. The relations to the HAs used throughout this paper 


are 


ttL/R — 4 / zjT/R 1 t jT/R\ 

H u±) = ^ H + ±H - ) 
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(C.ll) 


H S = H s , 


H t =H\ 


uTt — ^ ( rjTt I TjT t \ ttI t — zjTt 

■”||(±) — y 2 (•■"+ ^ > 




Hi = K 


In [35] the notation A^-, with i,j =||, _L, 0, is used for the transversity amplitudes. Note that 
when comparing to this paper the difference in the convention of the polarisation vectors has 
to be taken into account. 

C.4. G 1 ^ 1 for B —> K*Gl 2 in terms of Helicity Amplitudes for m^ 1 ^ m^ 

The formalism discussed in this paper allows a simple extension to the case ^ Tri£ 2 , so 
that the results presented in (C.10) can be adapted to test for possible lepton-flavour violating 
processes. Using the notation 

y2 = @1 Pi @2 @2 ’ 


(C.12) 


(A 7 * given in (B.6) and /3 P 2 = yfE p 2 ± 2 ), we obtain the following expressions for Gm’ le = 

J\TGm’ le (with M defined in (26)) : 


GS’° = \ {^E, + (Kl 2 + l#-T + |K? + (V A 

+ 4 m y ^ 2 (ki 2 + \hx\ 2 + ik i 2 ~(v^a 

+ 4 (eiE 2 - ^ |#T + ^ (eiE 2 + m ei m i2 + \H P \ 2 


H 


T t 


+ 


H 


T t 


+ 


H, 


16 ( A 7 * 

+ y ( 3 (E]E 2 + m tl mi 2 ) - 

+ l ( 3 {EiE 2 - m tl m t2 ) - (|K| 2 + \H T -\* + |K|' 


Tt 


KK + hYh 7 } + KK 


16 

+ y ( m,i 1 E 2 + mi 2 Ei) Im 

+ ^ {m h E 2 - m l2 E 1 ) Im + H*ff£ + KK] , 


&* = 


Re [KK - KK] + 2V2 m ^ ./^ Re 


+ 2 m. 1 +m la im 


K 2 




+ V 2 


Q 

m ei - mi 2 

vK 


H p H Pt - h t RE 


Im [KK - KK] 


-"^Re [H ( fn p ] - mtl + ™ 1 ' 2 Re [H%H s ] + Im V2H%H p + 2H Pt H s 




V? 


g°1 = (2 |K| 2 - |K| 2 - \ H - 1 2 + (y ->• a) - 2 (2 |K| 2 - |K| 2 - |K | 2 
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H Tt 

n o 


H 


T t 


H 


T t 


G 2 ’° = - 4 (W?: 


: 2 + hl 

9 V“ 4g 2 


Hl\ z + \HY\ Z - 2\H^\" + (V ^ A 


4m ^ (\HX I 2 + |Flf - 2 K| 2 - (1/ -> A 
3 


+ jj ( y E l E 2 - m h m h + yj) l^l' + !j (^ 1^2 + m tl m h + \H P \ 2 


16 / A 7 * 

~ y ( 3 (E 1 E 2 + m tl me 2 ) - ^ 




•T t 


+ 


H 


-T t 


- 2 


H, 


-T t 


(3(E l E 2 - m tl me 2 ) -h^)(\Hl\ 2 + \H P \ 2 -2 \ H?' 2 


16 

“ y (me 1 E 2 + me 2 Ei) Im 


4r 


- 2i7,] / //,? 


8\/2 


(vn^E 2 - m h E x ) Im + H -^- ~ 2H o^] , 


gI' 1 = - 


Re [#+#+ - + 2\Z2 m ^ 2 m ^ 2 Re 
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//I//l' - // 7 // T ' 


+ 2 (m^ +m <2 ) Im 


V 7 ? 


+ 2 


m R me 2 ^ _ r rrA - P 


- H^H 7 } 

me 1 + m£. 


+ mfa) Im [<i^ - 




L Re [H^H P ] +2- 


V 7 ? 


V 7 ? 

Re - 2 Im [V2^if p + 2 H Pt H s 


2 A. 


G'o 2 = (4 K| 2 + K| 2 + |^| 2 + {V->A)~ 2 (4 |tf 0 T | 2 + |if p | 2 + |FT| 2 ) 
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uT t 
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T t 
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-Tt 


G 2 { 1 = 4v ^y f (i^i# + H^hX - + mei t^ 2 ( HlH s + h s hY) 

v 3 \ vV 

- V 2 i (h p H p - H P H P + \pl [e s H t E - H P H‘ 

m h ~ me 2 (rjAuP , tjP frA \ 
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(H^H P + H P H f) 
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G 2 / = - 2 ( H'iH'i + 2 H^Hfj 


4 A 7 * 

3^ 

8 A, 


2 2 
— mi 

+ 2\/2———j—— 

r 

tv ttV , u-^fvy 


rr~i — rr~i /Ti — /t-t rr-i — rTi /Ti — rr~i 

+ BpB n T - H^H?} - 


^ + {V -> A) - 2 ( BjH,] + ^ ti 2 _ + 2 ( B^B 0 Tt + hJ;h t g 


tV fjV I [7 A frA 


t T tjT 


T T t uTt 


(C.13) 


C.5. Explicit Helicity Amplitudes in terms of Form Factors 

We collect here the definitions of the Helicity Amplitudes in terms of which our results are 
expressed. The hadronic HA is defined by 

Ff = (K*(\)\sT X b\B) , (C.14) 

—X 

with T' |a A '->-a as defined in Tab. 3.1 and the further replacement u -A Co . The definitions of the 
hadronic matrix elements used in the calculations are standard (e.g. [76]). Below we evaluate 
the HAs using form factors to make clear the relative signs between the various contributions, 
allowing for definite comparison with the literature. 

Results for form factors for low q 2 can be found from Light-Cone Sum Rules (LCSR) with 
vector distribution amplitudes (DA) in [5,76] and B-meson DA in [6], and for high q 2 from 
lattice QCD [77]. Long-distance effects contribute to only, and include quark loops (QL), 
the chromomagnetic operator Og, quark loop spectator scattering (QLSS) and weak annihila¬ 
tion (WA). At low q 2 , effects have been evaluated in QCD factorisation (QCDF) in the leading 
1 /rrq,-limit and in LCSR. Results for Og, WA and QLSS in QCDF are given in [7], and addi¬ 
tional contributions for Og in [8]. In Ref. [7] it was shown that quark loops can be integrated 
into the 1 /mf, framework using the results from inclusive matrix element computations [9]. 
Results for Og and WA, as well as a prescription for dealing with endpoint divergences of 
QLSS, can be found in [10] and [11]. Results for charm loops beyond the 1 /mt, approximation 
can be found in [12] for LCSR with B-meson DA, and [13,14] for LCSR (at q 2 = 0 only) for 
vector-meson DA. At high q 2 many of the long-distance contributions are suppressed in the 
formulation in terms of an OPE in 1 /q 2 (with q 2 ~ m 2 ) [78,79]. It should be added that the 
large contribution of broad charm resonances in B -A Kfi + /r - observed by the LHCb collabo¬ 
ration [80] demands a reassessment of duality violations [66]. Long-distance contributions can 
be found elsewhere. 

Explicit results for the B —> I\*(£ i^-inode are given by 


HX 

Ho A 

Hi 

Hi 


Aim. B m,K* [( Cy - C' v ) (m B + m K *) A 12 + m b (C 7 - C' 7 )T 2 3\ 
\fq 2 (m B + m K *) 


Aim B nriK* 

7 ¥ 


(C A - C'a)A\2 , 


o (m B + mj- ) + Cy)y/xl V - (■ m B + m K *) 2 ( C v — Cy)Alj 

+ ^ (±(CV + C' 7 )\f\ B T 1 - (C 7 - C' 7 ) (: m 2 B - m 2 K .) T 2 ) , 

2 ( mB ( ± <^ + C '^ v - ^ ~ c '^) ■ 
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{C A - C' A ^j A), 

s _iVX B (Cs^ + rne^m h ( ^ } ^ 

\m(, + m s <r / 


V 


mb + m s 


2\/2 mpmx* 
mp + mx* 
2 mpmx* 


(C r + C' r ) T 23 , 


^ (C T - C' T ) T 23 , 

mg + mp 

^L= (± (C r - C^) - (C'r + C' T ) (m| - <*) T 2 ) 


Hi* = 9 / 2 (± (Cr + C^) yXXTi - (Cr - Cf) (m| - <*) T 2 ) , (C.15) 


2yV 


where Cy( A ) = Cg( 10 ) in the standard notation used in the literature and the (^-dependence 
of the form factors is suppressed. Furthermore we have used 


A 12 = 

T'23 = 


(m B + mi-.) 2 (m| - rn 2 K , - q 2 ) Ai - X B A 2 
16 mpm 2 K „ [m B + mi*) 

(m| - ra^,) (m| + 3m|-» - q 2 ) T 2 - A B T 3 
8 mim|, {mp - mp) 


the same shorthand for zero-helicity form factor combinations as in [76,77]. 

The so-called timelike HAs, often denoted by Ht in the literature, have been absorbed 

into H s and H p . This is exceptional and follows from the vector and axial Ward identities 

q^u(l i) 7 /i [ 75 ]u(£ 2 ) = (me 1 =F mp 2 )u{i\) [75 }v(£ 2 ). A similar simplification procedure could be 

^— 

repeated by use of the equation of motion id p {sia^ v b) = — (m s + mb)s"/ f Jj + id^sb) — 2 si D M b 
(as used in [81]) for H J* if all of the operators present in the equation were used in the effective 
Hamiltonian. Since the higher derivative operators are not present in the effective Hamiltonian 
used in this paper, such a simplification does not occur. 


D. Specific Results for B —> K£ 


The angular distribution for this decay is 

d/^osde = G(0)jD °° 0 ( ^ } + G(1)jD °’° ^ + G(2)jD ° 2 0 m 


where, using the general leptonic HAs in appendix A.3 and taking me 1 7 ^ me 2 , the functions 
G^) = A(with M dehned in (26)), are given in terms of B —> K HAs by 


G(°) 


^4 (EiE 2 + me 1 me 2 ) + 



+ 


^4 (EiE 2 - m ei me 2 ) + 



+ 16 


(^E\E 2 + me 1 me 2 


A'y* 

12? 


h v | 2 + ^4 (EiE 2 - me L me 2 ) + 

) + (^(EiE 2 + m ll me 2 ) + 

j | h Tt | 2 + 8 ^E\E 2 - m fl m< 2 - 


h p 
I h T 


2 

2 
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+ 16 ( rri£ 1 E 2 + m^Ei) Im [/A /i Tt ] + 8\/2 ( m^E-i — rri( 2 Ei ) Im [/i A /i T ] , 


G (1) = -4^/V ^ 


m tl + 777.^2 ~ TOfe ^AfoP 


— Im 


2 h Tt h s + V2 h T h p 


vV 


G(2) = (J^| 2 + |^|2 _ 2 |^T|2 _ 4 IfcTt Q 


(D.l) 


The equivalent expressions for equal lepton masses are, using the notation G^ lt > = J\fq 2 G 
G (0) = \ (1 + 2 rh e 2 ) \h v \ 2 + ^ft 2 |/i A | 2 + 2/3| |/i s | 2 + 2 |/r p | 2 

O O 


+ - (l + &rh 2 ) \h Tt \ 2 + -f3 2 \ h T y + 16ru^Im [h v h Tt ~\ , 
3 3 


i 2 + \g 2 

G {1) = —4ft (2m* Re - Im 2/i Tt /i s + y/2h T h p ) 

G (2) = ^V|2 + |^A|2 _ 2 \h T \ 2 - 4 |h T f) , 

where we have used the shorthand rhi = m*/ . 


(D.2) 


D.l. Explicit B —»■ K Helicity Amplitudes in terms of Form Factors 

As for B —>• K*l\l 2 we quote the HAs for form factor contributions only, which allows for 
comparison with the literature. Form factor computations are available for low q 2 and high q 2 
from LCSR [15,16] and lattice QCD [17] respectively. Contributions to long-distance processes 
can be found in the same references as for the A'*-meson final state (quoted in appendix C.5). 
The form factor matrix elements relevant to B —> K transition, in standard parametrisation, 
are 


(K(p)\s'y lx b\B(p B )) = (pb +p)„ f+(q 2 ) + 


-9/u (/o(<7 2 ) - f+{q 2 )) , 


(K(p)\sa /lu b\B(p B )) =i {pb + p ]L Qu ~ ( Pb + p) u q^ 


fr(q 2 


J ms + mK 


{K(p)\sb\B(p B )) = 


m b - m s 


-fo(q 2 ) 


(D.3) 


with (K(p)\s'y /1 'f 5 b\B( pb)) = (K(p)\s'y 5 b\B(pB)) = 0 in QCD. The hadronic HA is defined by 

h x = (K\sT X b\B) , (D.4) 

—X 

where T |a x -70 as i R Tab. 3.1 with ui —> CD, containing the full set of dimension-six operators 
in the effective Hamiltonian (12). We find 


h v = 


\J A bk ( 2 m b 

2^/q 2 \mB + m K 


(C 7 + C 7 )f T + (Cv + C'y)U 
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(D.5) 


h A 

h s 

h p 


V A bk 

2^? 


(C A + C' A )f+, 


m 


B 


- m? K f (C s + C' s ) mg 1 - mg. 


mb - m s 


+ 


-( C v + C v ) 


m 


B 


-n± / (Cp + C» + + \ 

2 \ mb — m s q z J 


h T = —i 


h Tt = —i 


V^BK 


\/2 (■ m B + m K ) 


(C T - C’ T ) f T , 


V^BK 


2 (m B + m K ) 


(Cr + C^fT, 


fo , 
fo , 


where the Kallen function (cf. equation (B.l)) A bk = A(m^, rn 2 K . q 2 ) replaces A# = A (m 2 B ,m 2 K *, g 2 ) 
and Cy(A) = ^9(10) i R the standard notation used in the literature. 


D.2. Comparison with the literature 

The results for equal lepton masses (D.2) do agree with the results of reference [36] when 
6g —> it — 8g is taken into account. This is consistent with the angular conventions. In this 
paper we use the same conventions as LHCb [37], which differ from the ones of [36] by the 
transformation stated above. 


E. Ab —A (—> (p, n)ir) £i £ 2 Angular Distribution 

The decay — > A (— > (p, n)Tr) with a final-state proton or neutron, recently measured by 
the LHCb Collaboration [82], can also be considered within the generalised helicity formal¬ 
ism, and is particularly relevant because this decay can also be described using the effective 
Hamiltonian defined in A.2. In this case (5) becomes, in the rest frame of the A&, 


AI(Ha 6 , H*, HaIAa,,, A at, Ai, A2) ~ 

a 7 ,a Ai j 7 (E.l) 

= E (!!») Va(!!,) < AlAs , 

j J y 

where the leptonic HAs are the same as before and J\f\ N is the HA for the decay A —> Nn 
analogous to the (]k* Kn factor in the B —> K* decay, this time carrying non-trivial dependence 
on helicities owing to the final state particle N having spin- j,. The terms \ A are the HAs 
for the Ab —> A decay and can be again expressed in the form 

U XiXa = (A(A A )|sF V 6|A b (A A J) , (E.2) 

—X 

with the T the same as defined in Tab. 3.1. The resulting angular distribution can then be 
expressed as 


K(q 2 , Ha, fig) ~ Re 


r 0 , 0 o 0,0 


(H a , Q e ) + I # 1 !# 1 (H a , H £ ) + /C°’ 2 fi °’ 2 (H A , tog) 
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+ jcl’ 0 nl’° (n A , n £ ) + /c^ 1 (n A , n e ) + tc^n 1 / (n A , n e ) 
+ (n A , n*) + tc{’ 2 n {’ 2 (n A , n e )] , 


(E.3) 


where 11 A = (0,0 A ,O) and £l £ = ((f), 0 £ , —(f>). A theoretical angular analysis of this decay has 
been performed in [40,83]; in terms of the functions defined in [40], the ICm’ le above are 


K 


o’° = 3 (Kicc + 2 K laa ) 

1 


/Ci ’ 1 = K 3s + iIx As 


C 1 = K lc 


A, 


°’ 2 = - (K lcc - K laa ) 


Aq° = - (K 2cc + 2 K 2ssJ , 

-1,1 _ 


icY = k 2c , /cl’ 2 = ^ (k 2cc - K 2ss ) , 

K}{ 2 = (K-isc + iK isc ) . 


(E.4) 


These results can also be compared with those found in [24]; it follows that the MoM will be 
equally useful in future angular analyses of this decay. 
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